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We  developed  an  unique  experimental  system  to  perform  optical 
investigations  of  semiconductor  materials  in  environments  generally 
beyond  the  capability  of  other  techniques.  Using  the  solid  helium  pressure 
generating  method  and  a  superconducting  magnet,  we  performed  laser 
induced  luminescence  studies  on  a  variety  of  structures  at  low  temperature 
(<  4K).  high  hydrostatic  pressure,  and  high  magnetic  field.  A  quartz  fiber 
was  utilized  to  provide  an  optical  path  to  and  from  a  sample  subjected  to 
values  of  the  external  parameters.  The  use  of  the  fiber  alleviated  the  need 
for  pressure  chamber  windows  and  concern  for  the  transmissivitu  of  the 
pressure  medium,  and  allowed  investigations  over  a  frequency  range  limited 
only  by  the  optical  properties  of  the  fiber. 

Our  investigation  concentrated  on  the  effects  of  hydrostatic  pressure 
and  magnetic  field  on  strained  layer  superlattices,  specifically  on 

lnj(6a,_jjAs  /  GaAs  structures.  In  these  structures  the  value  of  x.  the 


fraction  of  In  present,  affects  the  amount  of  lattice  mismatch.  If  the 


layers  are  sufficiently  thin,  the  mismatch  is  accomodated  by  strain 
and  does  not  result  in  dislocations.  The  biaxial  components  of  the  strain 
remove  the  degeneracy  of  the  valence  band,  resulting,  in  effect,  tn  a  simple 
band  situation. 

This  interpretation  was  corroborated  by  our  magneto-optic  data  which 
showed  that  the  Landau  level  energies  followed  the  well  known  relation. 

and  yielded  a  value  for  the  hole  effective  mass  of  m^**  «  0.16mo.  This 


value  was  consistent  with  the  results  obtained  from  magneto-resistance 
measurements  (0.13  -  0.17mo)  and  indicated  a  conduction  band  to  light  hole 
band  transitioa  The  data  also  yielded  values  for  the  bandgap  energy  and. 
with  the  application  of  hydrostatic  pressure,  the  pressure  coefficient  was 
determined  to  be  *^12.2  meV/kbar,  in  good  agreement  with  the  values 
reported  for  m-v  materials  in  general  ('^12  *  2  meV/kbar).  The  pressure 
dependence  of  the  reduced  effective  mass  in  this  sample  was  found  to  be 


l.2X/kbar,  comparing  favorably  with  the  calculated  value  of  0.8X/kbar. 
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We  developed  an  unique  experimental  system  to  perform  optical 
investigations  of  semiconductor  materials  in  environments  generally 
beyond  the  capability  of  other  techniques.  Using  the  solid  helium  pressure 
generating  method  and  a  superconducting  magnet,  we  performed  laser 
^  induced  luminescence  studies  on  a  variety  of  structures  at  low  temperature 
(^4K).  high  hydrostatic  pressure,  and  high  magnetic  field.  A  quartz  fiber 
was  utilized  to  provide  an  optical  path  to  and  from  a  sample  subjected  to 
values  of  the  external  parameters.  The  use  of  the  fiber  alleviated  the  need 
for  pressure  chamber  windows  and  concern  for  the  transmissivity  of  the 
pressure  medium,  and  allowed  investigations  over  a  frequency  range  limited 
only  by  the  optical  properties  of  the  fiber. 

Our  investigation  concentrated  on  the  effects  of  hydrostatic  pressure 
and  magnetic  field  on  strained  layer  superlattices,  specifically  on 

In^Ga^.^As  /  GaAs  structures.  In  these  structures  the  value  of  x,  the 


fraction  of  in  present,  affects  the  amount  of  lattice  mismatch.  If  the 
layers  are  sufficiently  thin,  the  mismatch  is  accomodated  by  strain 
and  does  not  result  in  dislocations.  The  biaxial  components  of  the  strain 
remove  the  degeneracy  of  the  valence  band,  resulting,  in  effect,  in  a  simple 
band  situation. 

This  interpretation  was  corroborated  by  our  magnetO'Optic  data  which 
showed  that  the  Landau  level  energies  followed  the  well  known  relation, 

and  yielded  a  value  for  the  hole  effective  mass  of  m^**  %  O.I6mo.  This 

value  was  consistent  with  the  results  obtained  from  magneto-resistance 
measurements  (0.13  -  0.17mo)  and  indicated  a  conduction  band  to  light  hole 
band  ..  ansitioa  The  data  also  yielded  values  for  the  bandgap  energy  and. 
with  the  application  of  hydrostatic  pressure,  the  pressure  coefficient  was 
determined  to  be  *^12.2  meV/kbar,  in  good  agreement  with  the  values 
reported  for  III-V  materials  in  general  (*^12  ♦  2  meV/kbar).  The  pressure 
dependence  of  the  reduced  effective  mass  in  this  sample  was  found  to  be 
1.21^/kbar.  comparing  favorably  with  the  calculated  value  of  0.8X/kbar. 
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I.  INTROOIICTIOII 


Experiments  involving  the  variation  of  external  parameters,  e.g. 
pressure,  temperature,  magnetic  field,  have  long  been  used  to  determine 
optical  properties  of  solids.  We  have  developed  a  method  allowing  the 
simultaneous  low-temperature  optical  investigation  of  hydrostatic 
pressure  and  magnetic  field  effects  on  semiconductor  materials,  while 
avoiding  many  of  the  problems  usually  associated  with  experiments  of  this 
type.  By  attaching  the  sample  directly  to  an  optical  fiber,  we  are  able  to 
illuminate  the  sample  without  having  to  be  concerned  with  pressure 
chamber  and  Dewar  windows,  or  with  the  transmissivity  of  the  pressure 
medium.  The  luminescence  signal  from  the  sample  returns  through  the 
same  fiber  and  is  processed.  This  method  differs  from  more  standard 
pressure  experiments  where  hydrostatic  pressure  is  usually  applied  by 
immersing  the  sample  In  a  compressing  medium  -gas,  liquid,  or  soft  solid- 
and  optical  access  to  the  sample  is  through  this  medium  and  a  window.  The 
physical  properties  of  the  medium  and  window  limit  the  range  of  pressure, 
temperature,  and  frequency  at  which  optical  information  can  be  obtained. 
The  application  of  pressure  to  a  semiconductor  alters  the  lattice  spacing 
and  changes  the  energy  levels,  increasing  or  decreasing  the  energy  gap 


depending  on  material.  Since  the  electronic  properties  of  semiconductors 
depend  strongly  on  the  band  structure  and  the  value  of  the  forbidden  gap. 
marked  effects  are  expected  when  pressure  is  applied.  In  zincblende 
compounds,  for  example,  the  direct  gap  increases  with  pressure,  although 
GaAs.  beyond  a  certain  pressure,  becomes  an  indirect  gap  material. 

Subjecting  a  semiconductor  to  a  magnetic  field  causes  periodic  carrier 
motion  in  the  transverse  direction  and  results  in  quantized  orbits.  This 
raises  the  bottom  of  the  conduction  band  and  lowers  the  top  of  the  valence 
band  for  positive  effective  masses.  This  quantization  of  energy  in  the 
transverse  direction  changes  the  density  of  states  from  its  zero-field 
parabolic  configuration  to  discrete  (sub)  bands,  or  Landau  levels. 

The  semiconductor  samples  we  investigated  were  quantum  well  and 
super  lattice  structures.  The  term  super  lattice  is  applied  to  a  periodic 
heterostructure  of,  in  this  experiment.  Group  III-V  materials.  These 
structures  are  grown  by  a  process  known  as  molecular  beam  epitaxy,  which 
allows  accurate  control  of  layer  width  and  compositioa  Ttnis,  as  the 
number  of  layers  gets  large,  a  band-like  condition  analogous  to  the  Kronig- 
Penney  one-dimensional  crystal  is  approached.  This  impresses  a  periodic 
potential  on  the  structure  that  is  considerably  larger  than  atomic 


dimensions.  Since  the  bandgap  of  the  barrier  can  be  significantly  larger 
than  the  gap  of  the  well  material,  the  forbidden  gap  between  the  valence 
and  conduction  band  edges  will  exhibit  a  step  at  each  interface. 

For  barriers  of  sufficient  height  and  width,  an  electronic  carrier  can 
be  considered  as  confined  to  a  three-dimensional  potential  well.  If  the 
width  of  the  well  is  much  smaller  than  the  other  dimensions,  then  particle 
motion  in  that  direction  is  quantized,  as  opposed  to  the  unquantized  motion 
in  the  transverse  directions.  Thus,  in  the  heterostructure,  eigenvalues  for 
bound  states  can  be  found.  If  barrier  thicicness  however  is  reduced  to  the 
point  that  carriers  can  penetrate  the  barriers,  then  the  coupled  -  well 
superlattice  regime  is  reached.  When  the  wells  are  brought  together,  an 
interaction  between  the  degenerate  band  states  occurs,  causing  a  splitting 
of  the  the  well  levels.  This  splitting  is  a  function  of  barrier  height  and 
width. 

In  direct-gap  materials  band-to-band  transitions  between  states  of 
the  same  k-value  conserve  momentum  and  result  in  a  spectrum  given  by 
L(b)=K(ho  -  Eo)’^2.  where  K  is  a  constant  depending  on  the  reduced  effective 
mass.  Emission  therefore  has  a  low-energy  threshold  at  ho  =  Eo  and  a  high 
energy  tail,  in  addition  to  the  band-band  transitions,  electron-hole  pairs. 
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or  excitons,  may  also  contribute  to  the  emission  spectrum.  Excitonic 
spectra  are  characterized  by  an  ener^  below  Eq  and  a  narrow  linewldth. 

In  Section  H  we  develop  the  k*p  method  (bold  face  symbols  denote 
vectors,  except  «sh/2Tr)  and  expressions  for  the  energy  values  of  the 
bands  at  r  for  the  face  centered  cubic  lattice.  These  expressions  are  used 
in  Section  III  to  determine  the  direct  absorption  coefficient  in  the  <I00> 
directioa  We  also  discuss  the  matrix  elements  and  constants  and  address 
the  exciton  contribution  to  the  absorption  coefficient. 

In  Section  IV.  the  effect  of  hydrostatic  pressure  is  discussed  and  the 
absorption  coefficient  is  modified  to  include  the  pressure  dependence.  The 
magnetic  field  effects  are  discussed  in  Section  V  and  expressions  yielding 
the  Landau  levels  and  also  exciton  effects  are  presented. 

The  experimental  configuration  and  method  are  described  in  Section 
VI.  along  with  some  of  the  samples  used,  and  results  are  presented  in 
Section  VII.  Concluding  remarks  are  contained  in  Section  VIII. 
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II.  THE  k  -  p  HETHOD 

The  k  *  p  method  provides  a  powerful  tool  for  investigating  tr«  band 
structure  of  semiconductors  near  high  s;pmetry  points  in  k-space.  its 
advantage  lies  in  the  fact  that,  by  use  of  the  symmetry  properties  of  the 
crystal,  the  band  structure  at  a  point  k  depends  only  on  a  few  parameters 
which  may  be  determined  experimentally.  An  early  use  of  the  method  was 
to  find  an  expression  fo**  the  effective  mass  and  to  develop  the  "f-sum" 
rule.^'^  Subsequent  development  included  the  extension  to  degenerate 
bands^2)  and  consideration  of  spin-orbit  interaction.^^'^)  Although  the  k-p 
method  allows  calculation  of  band  structure  throughout  the  Brillouin  zone, 
it  is  most  readily  applied  at  selected  symmetry  points.  The  configuration 
of  the  empty  lattice  is  developed  and  shown  in  Appendix  A. 

In  order  to  develop  the  theory  of  the  k  •  p  method,  we  start  with  the 
one-electron  Hartree  Hamiltonian, 

H  ♦(  r)  =  l(-f|2/2m)  72  ♦  U(r)l  Hr)  =f  ’l’(r),  (I) 

where  U(r)  represents  a  periodic  potential.  The  eigenfunctions  of  the 
Schroedinger  equation  are  Bloch  functions  and  may  be  written  as 

♦rt(r)  =  c*'fUrt,(r),  (2) 


I  6 


where  n  represents  the  band  index,  tjf|((r)  has  the  periodicity  of  the 

potential  U(r).  andk  lies  in  the  first  Brillouin  zone. 

Substituting  E<^  (2)  into  Eg.  (f)  and  performing  the  indicated  operation, 
we  obtain 


l(-«V2m)v2  *  U(r)1e"‘f  UrtcCr)  =  ^Oc)«*Turt,(r). 


0) 


or. 


K,jUrtj(r)=l(p2/2m>(ft/mXk*pH*W2m)*U(r)lUrtt(r)  =  ^Oc)Urtc(rX  (4) 
Introducing  kg.  we  can  rewrite  Eg.  (4)  for  a  symmetry  point, 

I(p2/2m)  ♦  (f/m)(k  ♦ko  -ko)l)  ♦  (f2/2mXk2+ko2-ko2)  ♦  U(r)hJrti(r) 


and  separating. 


l(p2/2m)  ♦  (h/mXkoP)  ♦  (ft2ko2/2m)  ♦  Ur)  ♦  (h/mXk-ko)'  P 

*(h2/2mXk2  -  l(o2)lUrtfr)  = 


We  now  define  i  (p2/2m)  ♦  (h/mXkoP)  ♦  (f|2ko2/2m)  ♦ufr).  and. 
since  Mk^Urt^fr)  =€f^(ko)Urt5(r),  we  have 


^^0  *  *  (ftV2mXk^-ICo^)lUrt((r)  =  ^(k)Urtt(r).  (5) 

Writing  the  wavefunction  for  any  k  in  terms  of  ko. 

and  multiplying  both  sides  of  Eq.  (5)  by  and  then  integrating  over 

the  unit  cell,  we  get 

2nJsj^^rtto^''’'ko  *  *  (»>V2tnXk*  Hco*)!^'®'  •*«) 

«  €n(k)Ay^.(k  -ko).  (7) 

Performing  the  indicated  operations  results  in 

Sn'«n('^o)Snn'^¥^'('^**o)♦Zn•(^^^/2mXk2-ko2)8n„.ApJ^.(k-^ 

*Zn’ Jq  »/mXk  -ko>u"rt((r)pUrt(fr)^m•(k  ■''o)»  =  '•'o)  (8) 

Rearranging  and  dellning  Pnn’fto)  • 


we  obtain 


e 


♦  (ft/mXit  -ko)*Pnn'^o))Ann*Oc  **o) 


=€n^)Ajin'0c  -ko).  (9) 

Equation  (9)  is  the  eigenvalue  equation  in  the  ko  representation  for  the 
point  k  This  system  of  equatiorvs  can  be  solved  if  the  determinant  is  equal 
to  zero.  For  example,  in  the  case  of  a  simple  band,  index  n.  second  order 


perturbation  theory  yields 


*  (ft2/2mXk2  -ko2)  ♦  (A/mXk  -ko)Pnn*(ko) 

♦  (ft/m)22n*n‘^(‘^  -ko)*Pm'('^o)  (k  •*o>Pn'n(‘^oMfn<^o) "  (’0) 

As  a  further  example,  two  strongly  interacting  non-degenerate  bands  result 
in  the  matrix 


6p(ko)  ♦  (^iV2mXk2  -ko^)  -  e^Ck)  (f^/mXk  -koVPn'n^l^o) 


(f /mXk  -ko>Pnn'0co)  «n^o)  ♦(^>^/2mXk2-ko2)  - 


=  0.  (II) 


which  is  readily  solved  for  EpQc). 

The  extension  of  the  k  ■  p  method  to  degenerate  bands  is  developed 
below.  We  assume  an  s-fold  degeneracy  at  k  =  0.  The  set  of  functions 


a 


-S 

•  •  •  o 
1 


uio(0.  U2o(r).  U3o(rL..  UgoCr)  satisfying  HoUjo(r)  •  ^oUjo(r),  where  Ho  = 
(p2/2m)  *  iKr).  the  unperturbed  Hamiltonian,  belong  to  the  same  energy 


The  k  p  method  can  be  used  to  find  the  eigenfunctions.  Uji^Cr).  and 

eigenvalues,  £j(kl.  for  points  k^O.  As  a  zeroth  order  approximation 

we  let  Uji^fr)  =  A|(k)Ujo(r).  (12) 

Substituting  Eq.  (12)  into  Eq.  (5)  and  proceeding  as  before.  i.e.  multiplying 
by  u**j<o(r)  and  integrating  over  the  unit  cell,  we  have 

Ej  jQ^u*'j.o(r)M|jjjAj(k)Ujo(r)(l'  ♦  2j/Q^u*'j*o(r)(ft/mXk  -ko)'PAj(k)Ujo(r)dr 

•  Zj  Jq  u"cii(rX«V2tnXk2  -ko2)Aj(k)Ujo(r)(lr 

■ZjLu«j.o(r)£j(k)Aj(k)Ujo(r)(lr.  (13) 

Upon  separating  the  second  and  third  terms  and  integrating,  we  obtain  the 
following  expression 


Ij€j‘0co)Aj0c)8jj-  ♦Z](ft/m)(kpjj.  )AjOc)-2j(fi/mXlcoPjj'  )Aj(k)+ 

Zj(«V/2m)Aj(k)8jj.  -  Ij(*2koV2m)Aj(k)8jj.  =  6j.(k)A|(k).  (14) 

Rearranging  Eg.  (14), 

j:jl£j.(|co)$||.  -  (*2koV2m)8jj.  -  (*/mXko1>jj- )lAj(k) 

♦  Ijl(li2k2/2m)8jj.  ♦  (li/mXk  p|j.  )1A|0()  =£j.(k)AjOt).  (15) 

The  bracketed  I  ]  term  of  the  first  summation  in  Eq.  (15)  is  defined  as  €o- 
We  can  then  write 

{(^i/m)(k  pjj. )  -I6j.(k)  -  £o  -  (1s2k2/2m)]8jj.  )Aj(k)  =  0. 
or  simplifying  further, 

Z®j=,  t(«/m)(kDjj. )  -  EU)Sjj.  )Aj(k)  =  0.  (16) 

where  Etk)  »  Ej.(k)  -  ^  -(*W2m).  and  solutions  are  obtained  when  the 

determinant  of  the  system  of  equations  (order  s)  is  equal  to  zero. 

If  this  approximation  yields  zero,  it  becomes  necessary  to  go  to  the 
next  order  approximation.  This  is  accomplished  by  including  in  the  zeroth 
approximation  not  only  the  functior«  of  the  degenerate  band,  but  also 
others  at  k  »  0.  Using,  therefore, 


i|j(r)  =  Z®j=|  Aj(k)Ujo(r)  ♦  Z** |*jAj(k)Ujo(r),  (17) 
and  proceeding  as  before  we  get,  for  second  order  perturbation, 

[(f. Vm2)2 i(k  p jj  )(k  p j|.  )/(£o  j)  -  E"(k)S|j*  lAj.  (k)  =  0,  (18) 

where  E'’(k)  =  f 'Xk)  -  €o  -  (*2k2/2m),  and  f(k)  =  6o  ♦  f"(k).  Solutions  again 
are  obtained  when  the  determinant  of  the  system  of  equations  is  equal  to 
zero. 

As  an  illustrative  example  we  determine  the  perturbation  matrix 
for  a  crystal  of  cubic  symmetry  with  the  edges  of  the  cube  oriented 
along  the  x,  y,  z  axes.  We  have, 

-  «k)$jj.  lAj.(k)  =  0,  (19) 

where  0«>jj.  .(«2/2m)5|j.8„p*(*/m)2IiP«j|P>jj./(€o-«).  (20) 

with  <x,P  =  X.  y,  z  and  ^  the  origin  of  energy.  We  let  the  functions  *2. 
♦3  form  the  basis  of  representation  for  k  =  0,  which  transforms  under 
the  elements  of  the  point  group  of  the  cube  as  ♦,  ~  xfO').  *2  yf(ir),  and 
♦3  '^zf(r). 

Then,  s  (fi2/2m)8jj.  ♦  (l>/m)22jp^jjp*<jj./(£o-€i), 


and 


.  'f* . ' 


D^jj.  £  (li/m)2ZiP’‘]j  P^ij-  /(fo  -  fj). 

and  similarly  for  other  combinations  of  x,  y.  z.  The  matrix  elements  of 
perturbation  energy  are  thus  given  by 

■'♦j  I  *<kp  I  ♦  j’*  =  («'/m)^iWl  kl>  1  I  k-p  I  )’>V(  «o  -  <|).  (20 

Where  j,  j'  =  I.  2.  3.  and  i  ranges  over  all  other  bands. 

If  we  let  a:  a  constant,  then  from  completeness  of  the  set  we 
can  write 

*i  I  •‘kp  I  I  (kp)*  1 W  ^-«|>.  (22) 

Using  (kp)2=  li2lkk/8Vai<2)*k2y02/3t)k)*lt2202/822)*k,k||02/aiiay>...l  and 
r  =  («2  ♦  the  elements  of  the  matrix  mag  he  determined  one  by 


<1 1  (kp)2 1 1>  -  kJ Jxf(rX8V8x2Xxf(r))<k 

.  kk  Jxf(rX8k/8ykXxf(r)l(l'  •  „ .  (23) 

Upon  performing  the  operations  indicated,  we  find  that  the  coefficients  of 
kjkj  (I*))  are  odd  functions  and  thus  integrate  to  zero.  We  are  left  with 


Dll  =  ♦  (ft/m)2k2xl|l<l  |  P  x  1  I P  x  I  ^ 

♦(ft2/2m)(k2y  ♦k^p  ♦  (<i/mWy  j[<l  |  P  y  |  i><i  |  P  y  |  !>!/(  «o "  ^i)  • 

Expressions  of  similar  form  are  obtained  for  O22  ano  D33.  As 
examples  of  off-diagonal  elements  we  will  look  at  Dt2  and  D13. 

<1 1  (k  p)2 1 2>  ~  k2 Jxf(r )02/8x2)tyf(r)ldr 

♦  k2y/xf(r)OWXyf(r)lc^  ♦ ...  I. 

In  this  case  we  find  that  the  coefficients  of  k^j  integrate  to  zero,  as  well 
as  the  coefficients  of  kxk^  and  kyk2»  leaving 

0,2  =  (fi/m)2kxky2  j(<l  I  p  X I  i><i  |  p  y  1 2>l/(£o  “  €,) 

♦  (fi/m)2kxkyE  j[<l  I P  y  I  i><i  |  P  x  1 2>l/(€o  -  €,). 

For  0,3  the  result  is 

0,3  =  (li/m)2kxk2Sil<l  I  p  X I  |  P  zl  3>l/(^  -  €,) 

♦  (liWkxk^Z jl<l  I P  2I  i><i  I P  X 1 3>l/(€o  “  ^i)- 
The  corresponding  summations  in  the  above  two  expressions  are  equal  due 
to  symmetry.  The  remaining  elements  are  easily  found,  resulting  in  the 


S  k'k2*(c*(ft2k2/2m)  Bkyk^^lPf^kx  Bk^k^HP^^ky  Bk^ky+IPj^kj 


X 


D= 


Bkykz-iPnkx  LK^^+rKk^y+k^z)  Nlc^ky 


126 


Bk^kz-IPriky  NK^ky  LK^y+ttCk^x^'^^z^  ^Kykz 


♦  (ft2k2/2m) 


Bkx*^y"‘^n*^z 


N'k.k,  NK„k,  LKVM(k^«*k^ii) 


jji^Z  2  '  X  y' 


(^*(hW2m) 


The  coefficients  are  defined  as  follows: 


A'l  (ft/mFI|:'is|<S|p,|u|>|*/(fc- <l) 


B  •  2(ft/m)*2  jris<Slp,  |Uj><U|  |Pj,  |z>/l(l/2X  «c  •  «()  -  ^i' 
L-.  (*/m)2Z  j^i  I  <X  I  p,  I  u  j>  I  V(ey-£j) 

♦(«/m)^i'‘o  1  <X  I  p  J  ui>  I  V(  £v-£j)  =  F'  *26 


n  .  (*/m)2Z|*'iS|<x|py|u|>12/(£v-«l) 

*(«/m)*2|:‘25 1  <x  I  py  I  U|>  I  V(Cv-«j) 


N’  1  (f /m)2I  I  <x  I  p  J  U  j>  1 2/(£v-€j) 

-(l/2)(ft/m)22  I  I  Px  I 

♦(ft/m)I  j^B  I  <X  I  py  I  Uj>  I  V(£y-6j) 

-(<»/fn)2£  1^25 1  <x  1  Py  I  u  j>  I  V(«y-Cj). 

=  F‘  -  G  ♦  H,  -  H2 

In  each  case  the  summation  over  i  goes  over  higher  bands  of  symmetry 
type  shown  above  the  summation  sign.  The  conduction  band  (Ti  symmetry) 
and  the  valence  band  (Tfs  symmetry)  interact  strongly.  Weak  interactions 
with  other  bands,  according  to  Loewdin  perturbation  theory,  affecting 
coefficients  L'.  M,  N*.  lead  to  renormalized  interactions  between  the  bands 
of  the  strongly  interacting  set. 

The  matrix  of  Eg.  (26)  can  be  simplified.  If  we  limit  ourselves  to  very 
small  values  of  k  and  use  Loewdin  perturbation  theory  between  the  single 
conduction  band  and  the  three  valence  bands,  we  have,  considering  the 


interaction  Pf^  as  a  perturbation. 


vv^ 


S  Al(2  ♦ 


0=  X  0  Lk2^*  Nk^ky 


Nkxk2  .  (27) 


Nkjjky  Lk^y  ♦  M(k^j^  *  k^^) 

Nk^k2  Nkyk2  Lk22*lKk2x  ♦k^y)*  ^ 


where  the  coefficients  are  now  defined  as 


A  « (*/m)22  j^lS  1  <S  I  p  J  Ui>  I  V(€cq-€j)  ♦  PV(£cp-£^) 
L  s  I  <x  I  p  J  U j>  1 2/(€^-£j)  ♦  P2/(€^-€j) 

♦(fi/m)22  I  <x  I  Px  I U  j>  1 2/(€^-6,) 


n  1  as  before 


N«(l»/m)^  i^M  I  Px  I  I  ^A£yo"^i)*P^A€yo"^|) 

-(fi/m)2S  1^12 1  <X  |Px  I  Uj>  I  V(€yo-€|)  ♦  I  <X  I  Py  I  U|>  I  V(€y^-€|) 

I  <X|Py  |u,>|  2/(«^-f|). 

Until  this  point  the  effect  of  spin  has  not  been  considered.  The 


spIn-orbIt  Hamiltonian  Is, 


Kso-fi/4m2c2(7U>p)ff. 
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where  a  represents  the  Pauli  spin  operator.  Thus  the  total  Hamiltonian 
becomes 

^  total  *  ^-p  *  ^o. 

Including  Hgo.  as  given  by  Eq.  (28).  in  Eq.  (3). 

l(t»2/2m)72  ♦  iXr)  ♦  (^i/4m2c2XVU$')xp)-  a)i^(r)c’‘^'^  =^^j(r)e*'^'*’.  (29) 

or.  using  p  s  -iftO/ar). 
lpV2m  ♦  (^/m)(kp)  ♦  -h  V/2m  ♦U(r)li^(r) 

•»(h/4mV)(7U^')*p>cri^(r)-»(6V4m2c2XVU4-)«k)'  ffi^(r)  =  ^Ujj(r)  (30) 

Equation  (30)  can  be  rewritten  in  the  ko  representation  by  following  the 
steps  leading  to  Eq.  (5).  This  yields 
[pV2mHli/m)(kop)*fi2koV2m*U(r)Hli/m)(kHco)p*(fi2/2m)(k2-ko2) 

♦(ft/4m2c2)I7U^')><pl-  cr  ♦  (62/4m2c2X7Uj')*k)  •  ol  =  €n(k)Urt((r)  (31) 

or,  lH|jjj*(6/m)(k-ko)*  p*(62/2mXk2“ko2)K^»/4m2c2X7U^’)>«p)  •  o’ 

♦  (l•2/4m2c2X7U^’)>«k)•olUrtc(r)  =  ^Oc)Urk(r).  (32) 

where  ■  p2/2m  ♦  (ft/mXko’p)  ♦  hh^/2m  ♦Ur).  Defining 


Si 


u,^(r)  s  Zn'^nn'^~^o)Unicg^)'  ^  »Jt>stituting  this  e)q>ression  into  Eq. 

(32).  and,  as  before,  multiplying  both  sides  by  and  integrating  over 

the  unit  cell,  we  obtain 

*  («V2mXK’-ko’)15f„-  •  (*/mXlc-lco)  •  Pnn’ 


♦(«/4(nVXVlU’)«p>  a  •  («V4m2c»XVU«')>*)  ■  ®)Cnn.  =  «n*')Cnn’- 

The  addition  of  spin  results  in  eight  basis  states.  Hence  we  have  an 
eighth  order  matrix  and  therefore  an  eighth  order  secular  equation. 
Switching  to  a  new  set  of  basis  states  ,  the  spin  orbit  Hamiltonian  may 
be  written  in  diagonal  form.  The  states  which  diagonalize  the  spin-orbit 

Hamiltonian  are  those  corresponding  to  |j,  mj>  equal  to  |3/2.  3/2>. 

1 3/2,  l/2>.  and  |l/2.  I/2>(3), 

Although  two  additional  terms  were  introduced  in  the  wave  equation 
when  spin  was  considered,  only  the  k-independent  term  Is  kept  since  the 
crystal  momentum  is  very  much  smaller  than  the  electron  momentum.  The 
effect  of  spin-orbit  splitting  of  the  valence  band  (six-fold  degenerate 
p-states)  results  In  a  four-fold  degenerate  set,  Te  for  j  =  3/2,  and  a 
two-fold  degenerate  set.  Fj  for  j  =  1/2. 


The  eighth  order  secular  equation  may  be  solved  numerically  when  the 
values  or  the  unknown  parameters  are  taken  from  experiment.  In  order  to 
obtain  analytic  expressions,  certain  9>proximations  must  be  made.  For 
example,  the  k  *  p  interaction  matrix.  Eq.  (27).  is  accurate  for  small  kinetic 

energies  in  the  conduction  and  valence  bands  compared  to  ^  Adding  the 

spin-orbit  interaction  to  this  matrix,  we  obtain  an  matrix  composed  of 
a  2x2  block  corresponding  to  the  conduction  band,  a  4»4  block  which 
corresponds  to  the  valence  band  (j=3/2).  and  a  2  x  2  block  for  the  split-off 
band.  Diagonalizing,  one  obtains^^^ 

=  aV. 

for  the  conduction  band.  and.  with  Ag  representing  the  spin  orbit  splitting 
of  the  valence  band  and  given  by 

Ag  =  -(3ift/4m2c2)<XlOU(r)/a2)pj<  -  (9U(r)/ax)p2|Z>.  (34) 

for  the  valence  band 

Ey  3/2  =  €^0  ♦  As/3  ♦  ft2kV2m  ♦(L*2M)kV3 
t{(L-M)2k^4(1/3)lN2-(L-M)21(k2xk2y+k%k22*k2yk22)}’/2  (35) 


Ey  1/2  =  €vO  -  2V3  ♦  ♦(L‘*2n)k2/3. 


where  ♦(-)  in  E^  3/2  refers  to  the  light  (heavy)  hole  band,  and 

A-sA**Pfi2l2/3Eg*  l/3(Eg*As)l 
L-  .  F*  ♦  26  -  P„2/(Eg*As). 

with  Eg  1  €^-(€yo  *  ^  constants  L.  M.  N.  A*.  P|^.  F*.  and  6  as 

defined  earlier. 

In  the  narrow-gap  approximation,  applicable  to  GaAs,  the  denominator 
of  the  constants  6,  and  H2  is  much  larger  than  the  denominator  of  F,  and 

the  s<iuare  root  in  the  expression  for  Ey  3/2  can  be  expanded.  Rewriting  and 
neglecting  terms  of  second  order,  we  get, 
Ehh*Eyo*As/3*fi2k2/2m+Mk2HL-ri-NXl(2  xk2y*k2^k22*k2yk22)/k2  (36) 

E,hSfEyo*Ag/3+f2|f2/2m*(2L*n)k2/3-(L-n-N)(k2 
where  the  heavy  hole  band  corresponds  to  the  quantum  number  mj=i3/2,  and 
the  light  hole  band  to  mj=il/2. 

In  this  approximation  where  k  was  assumed  to  be  small,  the 
conduction  and  valence  bands  were  treated  according  to  Loewdin 


perturbation  theory  for  the  k  -  p  interaction,  and  with  spin  introduced,  we 


obtain,  for  a  small  direct-gap  material,  the  following  expressions  for  the 
energies: 

Ec«£cO  ^  V/2m^A*k2*p„2k2l2/3Eg^l/3(E  g+Ag)! 
Ehh*^*^s/3**^^/2m*nk2*(L-n-N)(k2  xk^y^k^^k^^+k^yk^^Vk^  (37) 

E,h«€^^As/3+li2k2/2m4(2L^n)k2/3-(L-n-N)(k%k2y+k2^k22+k2yk22)/k2 
Eso*^v0-2As/3+^i2k2/2m*(L*2t1)k2/3+P  g-1/3(E  g+Ag)!. 

Since  €yo  represents  the  k  =  0  valence  band  energy  without  spin-orbit 
energy,  we  replace  C^o  ♦  Ag/3  by  Ey,  and  €yo‘2As/3  by  Ey-Ag. 
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III.  THE  DIRECT  ABSORPTION  COEFFICIENT 
A.  <I00>  Direction 

We  considered  the  case  of  absorption  in  the  <100>  direction  and  also 
calculated  the  coefficient  for  the  directionally  averaged  case.  We  used  Eqs. 

(37)  to  calculate  the  joint  density  of  states  with  corresponding  to  E^ 

and  E^^.  E|^.  E^^  corresponding  to  Ej  in  Eq.  (B-21).  We  first  examined  the 
heavy  hole  to  conduction  band  transition.  We  have 

Jc.hh  =  JB2«k(2/8n3)8{Eg-li<o*A*kM(2PjiV3Eg)*P,^V3(Eg*As)llc2-mc2 

-(L-M-N)(kx2ky2^k^2k22*ky2|(/)/k2}. 
or,  in  spherical  coordinates, 

■’c.hh  =  JBz(2/8Tt’)8feg-t(i>HA'«(2Pf,2/3Eg)*P„2/3(Eg*As)-”*<* 

-(L-ri-N)k2(sin‘^ecos2esin2t+sin2ecos2e))k2sin0dkded».  (38) 

In  the  <100>  direction  0  =  (7r/2),  ♦  =  0,  and  Eq.  (38)  becomes 

Jc,hh  =  (2/87r3)Jg2S[A  ♦  Bk^lk^dkded*.  (39) 

where  we  defined  A  s  E«-li<i) 


B  .  A*H2P„V3Eg)*P„2/3(Eg*As)“M. 

Letting  X  =  A^Bk^,  then  dX  =  2kBdk.  and  integrating  over  e  and  we  have 

Jc,hh  =  (2/8it3)(4tt)  J  8(X)k2dX/21(B. 
or,  substituting  for  k  in  terms  of  X. 

Jc.hh  =  (2/87r3)(47T)(l/2)Js(X)(X-A)’/2/B3/2 

=  (27T2r‘B‘3/2(-A)»^2  (40) 

This  result  was  substituted  into  Eg.  (B-20).  and  yielded,  for  the  <100> 
direction, 

o<hh-c(“M2e2/ncin2u)  |  21A'*2Pf,2/3Eg*Pf,2/3(Eg*As)-Mr’'2 

« (»<1>-Eg)’^2  (41) 

The  same  method,  applied  to  the  light  hole-conduction  band  transition, 
resulted  in 

<>c,n.c(<o)a(2e2/ncm2<a)  |  K’,n-c  I  ^A'*2Pf,2/3Eg*Pf|2/3(Eg*As) 

-(2L*t1)/3r’'*«  (Aia-Egy/^,  (42) 

and.  for  the  split-off  to  conduction  band  transition. 


o<so-<J(<«W2e^/nctt|2m)|H'so.cl*lA'‘2Pf,2/3Eg-<L»2l1)/jr’^2 

-(ftia-Eg-V^*  (43) 

Combining  Eqs.  (41)-(43)  to  account  for  the  summation  in  Eq.  (B-20),  we 
have 

<100>  :  <x(<0)  =  *  «SO-C(6>) 

When  the  directional  dependence,  as  appearing  in  Eq.  (38).  is  isotropically 
averaged,  the  coefficient  of  the  anisotropy  term,  (L-M-N),  becomes  0.2,  and 
the  components  of  the  averaged  direct  absorption  coefficient  are 

«f(,.c(6)) » (2B2/ncm2«)  I  H^-cl  2|A-*2P„»/3Eg*Pf,2/3(Eg*As) 

-I1-0.2(L-(1-N)r’^^  » (*<i)-Eg)'^2  (44) 

«ir,-c(“)  “  (2e2/ncm2u)  |  K|h-c  I  *  (A’*2P„2/3Eg*P„2/3(Eg*is) 

-(a‘M)/3-0.2(L-n-N)P^2  X  (*<D-Eg)'/J.  (45) 

The  contribution  of  the  split-off  band,  being  independent  of  direction,  is 

still  given  by  Eq.  (43). 

B.  flatrix  Elements  and  Constants 

Examination  of  the  expressions  for  o((a>)  shows  that  they  contain  the 


following  constants  which  were  defined  in  Section  II:  A',  L,  tl.  N,  P|^.  They 
may  be  expressed  in  terms  of  the  parameters  6p.  Cq.  for  diamond 
structure  materials,  where  (p 

fQ'  I<rs|p|r2s>l* 

«R'’  l<r-B|p|r'2s>l*. 

with  r'25  representing  the  valence  band.  In  terms  of  these  parameters  the 
numerical  values  of  the  matrix  elements  for  Group  IV  materials  are. 

F  =  -fp/Eg 

ii=-<Q/i{(rB)-£(r2s)i  (46) 

6  =  "€p/l6(r’i2)"t(r'25)l 
For  Group  IV  materials  the  values  of  Cp.  Cq.  €p  are: 

£p  =  25.9  ev.  €q  =  17.2  eV.  €r  =  14.2  eV. 

Equations  (46)  are  transformed  for  the  zincblende  structure  by  including  an 
anti-symmetric  potential.^®^  This  results  in  Eqs.  (47). 

F  ="£pl£(n5^)-6(r^5y)'»£(r|5)-€(r’25)]/2U(r) 


M  =  -fQ/ff  (r  I5c)"f (I’  15 

G  =-6pi£(ri5j;)“6(rtsy)*£(r  is)*£(r*25)l/2i€(r  ^2  c)"^(rt5y)ii£(ri5c)-£(ri5y)i. 


and,  for  the  Interaction  between  the  conduction  band  and  the  higher 
state. 

A'  =  ”£Pl£(rjs^)-£(r|5y)-£(rts)‘»£(r*2s)l/2l£(r|5  ^)-£(r^5y)-Egl 

U(rt5cH(ri5v)l  (48) 

In  Eqs.  (47)  and  (48).  r  with  subscripts  c.  v  refers  to  the  Group  Ill-V 
material,  while  r  without  the  letter  subscripts  refers  to  the  corresponding 
Group  /V  material.  If  the  energy  levels  at  r  are  Icnown.  e.g.  from 
experiment,  the  constants  can  be  calculated.  Using  the  representative 
values(®^  for  Ge  and  GaAs  ( in  cV ), 

Ge:  £(r*2s)  =  0.  £(ri5) »  2.8,  £(r'i2)  *  12 
GaAs:  £(r^5y)  *  0,  £(rt^)  *  1.5,  £(ri5^)  ®  4.6,  c(ri2c)  *  12, 

then  the  values  of  the  constants  are 

F  a  -13.89.  M  a  -3.74.  G  a  -0.95.  L  a  -15.79.  N  a  -16.68,  A*  a  -1.63. 
These  results  can  can  be  checked  quickly  by  using  Eqs.  (37)  and  calculating 
the  effective  masses.  For  the  directionally  averaged  case  (T->0)  we  have 


2/£.l£/(U20)-(®V*®3)£/2‘^d2^Vfe  1  (0)2UJOU/2a2)«(c«)]p-®s>o 


W  miM 

ue6aq  aM  ‘uoipajip  jo  ssaipjeSaj  ujjoj  auies  au)  seq  pueq  ijo-mds  aiu 
JO  uojjnqijjuoo  am  aouig  ’paijiiduijs  ajaM  suoijoajip  snoiJCA  aqj  uj 
(co)>o  joj  suojssajdxa  au)  ‘sjuejsuoo  ag)  jo  sanieA  pajeinoieo  agj  6ujsn 

anieA  leoiJauinN  o)  («)»  jo  uo||onpau  *3 

‘papaiDau  aje  sngj  pue  ‘uouejedas  jaSjei  oi 
anp  jaiieujs  gonui  jaAaMog  aje  Ragi  ’jsjxa  spueq  jag6|g  uiojj  suoijnqijjuoo 
gOnogge  ‘(g^)  pue  sba  m  ubaiD  se  spueq  uaaMjaq  suogoejajuj 
asogj  0)  saAiasjno  paguiti  nueuouuaju!  pm  sjuejsuoo  agj  BugenieAa 
ui  juapjAa  si  juauiaajOe  aiqeuoseay  '^g^sjuaujijadxa  snoijeA  ujojj  paAjjap 
sanieA  pajdaooe  fijuoujujoo  aje  sjsagjuajed  uj  uaAiS  sanieA  uosjjeduioo 

'(®uiwro)  °uJ93ro  « 

(OuJ380'0)  °UJS80’0« 

(0U1|9  0  )  ®UJIQ9'0  as  ‘^«UJ 
(Ouj^900)  ®UJ9900ss 


where  H'  =  <f  |c*p  1 1>  =  c-<f  |  p  1 1>  =  cose<f  |p  |  i>.  9  represents  the  angle 
between  the  polarization  vector  and  the  dipole  moment.  Since  we  are 

concerned  with  transitions  between  the  (valence)  band  and  the  Fi^ 
(conduction)  band,  we  can  relate  the  matrix  element  of  K‘  to  P|^.  and  hence 
to  the  parameter  fp,  |  H'  p  =  cos^e  |  <Fi(.  |  P  |  Fts  y>  P 

=  cos^e  (49) 

=  cos20(m€p/2). 

Using  the  calculated  values  of  the  matrix  eiements  appearing  in  Eq.  (43) 
above,  the  mass  term  can  be  calculated. 

lA'*2Po2/3(Eg*A)-(L*2n)/3J’5/2  =  (I5.43)'5^2  =  o.oi65.  (50) 

Substituting  Eqs.  (49)  and  (50)  into  Eq.  (43),  and  remembering  that  the 

expression  (50)  is  in  terms  of  *V2m.  <x5q_^X(i>)  becomes 

«5o-<^^)*0.0165(2^  ^Vcp/  m/ncfi^)cos20u>"’(fio>-Eg-As)'^*. 

Checking  dimensionality  we  get 

«so-c?^)  KE^D^E*/ E/cl/lcx(E2>‘t2)x/ El  -»-» L‘\  as  required. 

Then, 


«so.^X<dM0.0165)  {123/2(4.8x10  "*0)225.9(1.6x10  "*2X9.1x10  "2®)*/2 


(l.3xl0"0)/t3.4(3xl0*0Xl.06xl0-22)23(1.6xl0-*2)l  K*(D)"’(li<D-Eg-As)*/2 

where  the  average  value  of  cos2e  was  taken  as  1/3  and  n,  the  refractive 
index,  as  3.4.  Performing  the  calculation  yielded 

«SO-<^0>>«(00>65X2.I63x|0SX  6(D)"*(f(D-Eg-As)*/2  (51) 

Similarly  we  can  write  and  o(|^_^(<i>)  for  the  direction  considered: 

<100>:  e<j^.c(0))s<0.0130X2.163xl0«X  60))"*(l»a)-Eg-As)*/2 

«lh-c(0)>8(0.0075X2.163xi0«X  6(o)"*(fi(i)-Eg-As)*/2  (52) 

<Avg>:  (xj^.c(<i))K0.0140X2.163xl00X  fia>)"*(fiCD-Eg-As)*/2 

e<|h-c(<«>)K0.0071)(2.163xl06X  6<i>)"’(6a>-Eg-As)*/2  (53) 

Appropriately  combining  Eqs  (51)  -  (53)  gives  the  expressions  for  the 
absorption  coefficient  for  the  selected  directions  in  the  crystal.  This  was 
done  in  Eq.  (54). 

<100>:  o<(<»))«44342(6<o)"*(<^(«)-Eg)*/2+35960(6a>)"*(6a)-Eg-As)’/2 


<Avg>:  «((o)*45639(liU))"*(fia)-En)*/2*35690(fi(«>)"*(li(i>-En-Ac)*/2.  (54) 
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•*  *1*1'*^^  '*iT‘*  i'^  ^ 


Units  of  o((<i))  are  cm~'. 

D.  Eiciton  Contribution  to  «(«») 

Until  this  point  the  effect  of  excitons  has  not  been  considered  in 
developing  expressions  for  the  absorption  coefficient,  hence  designated 

(xo((o).  o(o((i>)  was  shown  to  equal  zero  for  values  of  f  (i><Eg,  and  increased 

as  (t%(o)"Kt»6>-Eg)^'^2  for  Experimental  evidence  however  suggests  a 

more  complicated  relationship  between  the  absorption  coefficient  and  the 
incident  light.  Specifically,  when,  in  the  absorption  process,  an  exciton  is 
created.^^®'^’^  the  Coulomb  interaction  between  the  electron  and  hole  can 
significantly  affect  the  absorption  coefficient,  especially  near  the 
absorption  edge.  Applying  Elliott's  theory.(i2^as  developed  for  simple 
non-degenerate  parabolic  bands,  to  the  more  complicated  case  of  GaAs.  we 

find  that  the  absorption  coefficient  for  t^(i>>Eg  is  modified  by  a  factor 
27r2fj{/[1-exp(-27r3j,)l.^’^‘’^^  Thus,  in  the  region  of  continuous  absorption. 


«(<!))  *  «o(<d)  *  2Tr3„/[l-exp(-27r3  „)!, 


where  3^  «  (Rxi/(fta)-Eg)l’^2  and  the  value  of  the  exciton  binding  energy  In 


GaAs  Is  reported  to  be  In  the  range  R^i «  4.2  -4.4  mevJ 
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Equation  (55)  clearly  shows  that  the  exciton  effect  near  fio)  «  Eg  is 

significant,  diminishing  as  increases.  When  or  4^<i)-Eg»Rjj|. 
then  o(((i>)  ->  o(o(6>).  as  can  be  seen  by  expanding  the  denominator  in  Eq.  (55) 
and  neglecting  terms  of  order  higher  than  Conversely,  for  ^ 

fiti)-Eg«Rxi.  then  «((!))  -» o<o(o>)2ir2rj{,  where  the  exponential  factor  is 
considered  to  be  approximately  zero. 

The  quantity  R^i  is  related  to  the  gap  energy.  Eg,  and  the  exciton 

energy.  E^.  by  the  relation 

Ej,  =  Eg  -  Rxi/n2 .  n  =  1.  2. 3 . 

and  Elliott^’^)  predicted  a  series  of  lines  at  energies  lia>  =  Eg  -  R^/n^,  with 

intensities  decreasing  as  n-^.  While  for  GaAs  the  first  of  these  lines, 
corresponding  to  the  lowest  exciton  state,  occurs  at  an  energy  '^4.3  meV 

below  Eg,  subsequent  lines  are  much  closer  and  can  be  considered  to 
consist  of  a  quasi-continuum  with  oc((i>)  «  27r2fj^<xo(a>).  thus  transitioning 
smoothly  to  o<(cd)  for  ll(»)>Eg, 


Since  the  lines  for  n>2  are  not  usually  discernable.  the  absorption 
spectrum  at  low  temperatures  generally  shows  a  sharp  peak  at  -  Eg  -  R^|. 

for  GaAs  «(E^)  ^  2.5x10^  cm"^  and  decreases  to  the  approximately 

constant  value  near  Eg  before  approaching  o<o(o>)  for  fio)»Eg.  Thus,  in  Ill-V 

compounds,  the  primary  effect  of  the  exciton  on  the  absorption  coefficient 
is  to  modify  its  value  near  the  band  edge. 


IV.  HYDROSTATIC  PRESSURE  EFFECTS 

The  application  of  hydrostatic  pressure  to  a  crystal  causes  a  change  in 
volume,  i.e.  it  t)rings  the  atoms  closer  together.  The  decreased  separation 
between  atoms  results  in  a  changed  energy  gap  and  thus  the  separation 

between  two  bands,  changes  with  pressure.  Hydrostatic  pressure 

does  not,  however,  change  the  crystal  symmetry,  and  hence  the  selection 
rules  for  the  transitions  remain  unchanged.  In  addition  to  the  constant 

matrix  elements,  due  to  spin-splitting  may  also  be  considered  constant 

with  pressure.^'^)  A  small  change  in  the  lattice  constant  causes  a  linear 

variation  in  E^,  and  hence  we  can  write 

En  =  EpCO)  ♦  (dEn/dP)TAP  (56) 

where  (dEf/dP)j  denotes  the  pressure  coefficient  at  constant  temperature. 
The  coefficient  is  different  for  the  various  bands  and  may  be  positive  or 
negative.  (In  zincblende  structures  dEg/dP  >  0,  i.e.  Tg  moves  up  with 

respect  to  Tg).  Values  of  the  pressure  coefficients  for  Group  IV  and  Group 
lll-v  materials  change  only  slightly  from  one  substance  to  another,  and 
accepted  values  for  the  direct  gap  (zone  center)  for  GaAs  are:^’^“22) 


(57) 


(dEo/dP)x  »  GV/bar 

(dE,/dP)T  s!  S^IO"®  GV/bar 
(dE2/dP)T  «  3.6*10"®  GV/bar, 

whGTG  Ej  (i  =  0, 1,  2)  rGfGrs  to  bands  r,5^  and  r,2j.  rGspGCtivGly. 

Having  GarliGr  obtainGd  numGrical  valuGS  for  thG  constants  F.  M.  G. 
Gtc.,  WG  now  modify  thGm  to  includG  thG  GffGCt  of  hydrostatic  prGSSurG. 
SincG  thG  matrix  GiGmGnts  thGmsGivGS  arG  considGrGd  invariant  with 
prGSSUTG,  WG  havG,  for  GxampiG. 

|<  I  >1  F£o  =  Fp(Eo+AEo pX 

WhGTG  thG  subscript  P  dGnotGS  thG  prGssurG  dGpGndGncG,  and  from  Eq.  (57) 
wGhavG 

p  =  OEo/aP)jAP. 

Similarly  wg  obtain  for  the  other  "constants', 

rip  =  nE|/(E,*AE,  p). 

Op  =  GE2/(E2'^AE2  pX 

A'p=A‘(E^“Eo)/l(Ei*AEt  p)-(Eo'»AEo  p)l. 


Upon  substituting  the  new  constants  into  the  expressions  for  o(((i>).  we 


obtain  o(((i>,P).  For  example,  we  can  write  the  split-off  band  contribution 
to  the  pressure  dependent  absorption  coefficient  as 


i 

i 

«so-<^^‘P)*2.I63>‘106(A‘  p*2P2„/3(Eo+AEop*AsHLp*2np)/3l“3/2 

H(*<i))‘’(f»<D-As-Eo-AEop)V2  (58) 

where  P^^  continues  to  denote  the  matrix  element  and  P  the  pressure. 

Equation  (58)  can  be  evaluated  readily  if  pressure-modified  constants  are 
calculated  and  used,  or  it  may  be  recast  into  a  form  including  the  pressure 
coefficients  directly.  Thus, 

«so.<J[(i)f)«2.163x|0‘  {A'(E,-E(,)/HE|*M,_p)-(Eo*AEo,p)l 
*2P2|i/3(E««AEo_p*As)-FE<|/3(E(i*AEo_p)-2eE2/3(E2«AEjp) 

-2ME/3(E,»AE,P)}-A*4>rK*U-As-Eo-AEo_p)''2,  (59) 

Where  the  constants  A*.  F,  G,  and  n  have  the  values  listed  after  Eg.  (48). 

Although  relevant  information  Is  contained  in  Eq.  (59),  It  may  be  presented 
in  a  more  appealing  form  after  some  algebraic  manipulation.  Expanding 
applicable  terms  of  Eq.  (59)  and  limiting  ourselves  to  first  order,  the  terms 
are  modified  as  follows; 


A'(E|-Eo)/ (E  rEo*E  1  p  “Eo,p ) 

2PV3(Eo*AEop*As) 


AH'C  AE|  p  -AEq  p)/ (EpEo)] 


-  l2P2„/3(Eo^As)][I-AEo  pAEo^Ag)! 

F(o/(Eo^AEop)  -  F[hAEop/Eol 

2GE2/3(E2>AE2p)  -  (2G/3Xl-AE2p/E2] 

2ME,/3(E,*AE,p)  -  (2n/3)ri-AE,p/E,I. 

Substituting  these  approximations  into  Eg.  (59)  yields. 

«so.^J(6>.P)^2.163x|06  {a11-(AE,  ^-AEq  p)/(ErEo)] 
A2P^„/3(Eo^As)ffl-AEo  pAEo^Ag)?  -  (F/3)fl-AEo,p/Eo] 
-(2G/3)fl-E,p/E,H2ri/3)ri-E,  ,p/E,l}‘3/2 

Alia>)"K*a>-Ag-Eo-AEo  p)’^^. 

Separating  the  terms  in  brackets  (  }  into  pressure-dependent  and 
pressure-independent  groups  and  defining 

Q  ■  A'*2P2,^/3(Eo*Ag)'F/3-2G/3-2M/3 


R(P)  >  -AT(AE,p-AEop)AErEo)]-2P2„AEop/3(Eo*Ag)2 
♦FAEo  p/3Eo*2GAE2. p/3E2*2riAEi  p/3E,. 


(60) 
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we  can  write 

0(so-<^<^.P)*2.163xl06lQ*R(P)l  '^2(ft(o)'Kli<D-As-Eo-AEop)’^2.  (51) 

Again  expanding  to  first  order  we  have. 

«so-c^^‘P^*2.163x10SQ  ■3''2ji_3r(p)/2q](  fia>rKi»a>-As-Eo-AEop)^''2  (52) 

It  is  possible  to  write  Eg.  (62)  in  terms  of  O(o(o>)  by  simply  expanding 
the  square-root  term  and  then  combining  the  two  truncated  series,  keeping 

only  terms  up  to  first  order  In  AEj^  The  result  is 

■3/2(fa))-»(fi<o-As-E)’/^[l-3R(P)/2Q 
-  AEo p/2(li(D-As*E)*  h.o.t.)  (63) 

which  has  the  form 

«so-c^<^.P)  *  0(o(<i))>‘f(P). 

The  quantity  E  in  Eq.  (63)  in  this  case  must  be  taken  as  Eq  ♦  AEop  in  order 
to  account  for  the  pressure-induced  shift  in  the  energy  gap.  It  is  apparent 
that  Eq.  (63)  reduces  to  the  expression  obtained  earlier  for  ®-  9’ 

Eq.  (51).  when  AP  =  0.  The  contribution  of  the  split-off  band  to  the 
pressure  -  dependent  absorption  coefficient  can  thus  be  calculated  from 


Eqs.  (58)  or  (59),  or  using  the  Definitions  (60),  by  Eqs.  (61)  -  (63).  We 
prefer  the  form  of  Eq.  (62)  and  will  use  this  form  for  the  remaining 
contributions. 

The  contribution  of  the  heavy  hole  band  in  the  <100>  direction  can  be 
written  as 


«hh-c(^-P)  *  2.163x106{A‘l1-(AE,  p-Mop)/(ErEo)l 

*(2P2„/3Eo)ll-AEo  p/EoWP2„/3(Eo*As)l[1-AEo  pAEo^Ag) 
-Ml1-AE,p/E,ir3%6>)'Kli6>-Eo-AEop)'''2,  (54) 

where  the  expansions  are  again  to  first  order. 

By  the  same  procedure  as  used  in  the  case  of  the  split-off  band,  Eq. 
(64)  can  be  put  into  a  more  convenient  form.  The  result  is 

«hh-c(<*>'P)*2.163x|08Q  ■3/2[i-3r(p)/2Q1(  li(i))-»(fi(i>-Eo-AEop)’''2  (65) 

where  now 


Q  »  A’-2P2„/3Eo*P2„/3(Eo^As)-M 
R(P)  s  -A*(AEip-AEop)/(Ei-Eo)-2P2„AEop/3Eo2-P2nAEo.P/3(Eo*A  5)2 


I 


For  the  light  hole  band  we  obtain 

o<,h-c(o>.P)»2.163>‘10SQ  ‘5''2t^-3R(p)/2Qi(  fta))"Ki^(0-Eo-AEo)’^2, 

which  is  identical  to  Eq.  (65).  but  where  the  definitions  of  Q  and  R(P)  are, 
for  the  <100>  dlrectioa 

Q  s  A‘-»2P2„/3Eo*P2m/3(Eo*As)'2F/3-4G/3-I1/3 
R(P)  =  -A‘(M,  p-AEo  p)/(ErEo)-2P2„AEop/3Eo2-p2„AEop/3(Eo*As)2 

♦2FAEop/3EoMGM2p/3E2*nAE,  p/3E,  (67) 

Within  our  approximations,  the  general  expression  for  the  absorption 
coefficient  will  be  of  the  form  of  Eq.  (65).  with  0  and  R(P)  varying  as  we 
consider  different  directions.  For  the  <Avg>  direction,  heavy  hole 
contribution,  we  use 

Q  =  A*+2P2,^/3Eo*P2f^/3(Eo*As)-0.60G-0.60ri 

R(P)  s  -A*(AE,  p-AEop)/ErEo)-2P2„AEop/3Eo2-p2„AEop/3(Eo*As) 

♦0.60GAE2p/E2*0.60riAE,  p/E, .  (68) 

while  for  the  <Avg>  direction,  light  hole  contribution,  we  use 


Q  s  A*-*2P2„/3Eo+P2fi/3(Eo*As)*2F/3-0.73G-0.73n 


R(P)  s  -A*(AE,  p-AEo  p)/(ErEo)-2P2„AEo  p/3Eo2-p2„AEo^/3(Eo*As)2 


♦2FAEo  p/3Eo*0.73GAE2  p/E2*0.73nAE,  p/E,  (69) 

in  order  to  obtain  o<(c;>.P),  the  directionally  independent  split-off  band  must 
be  added  to  the  appropriate  «hf^_^a),P)  and 

Since  hydrostatic  compression  decreases  the  lattice  constant  and 
increases  the  width  of  the  (direct)  energy  gap.  experiments  have  been  done 

to  determine  the  pressure  dependence  of  Eg.  Welber  and  coworkers^^o) 

applied  pressure  to  GaAs  at  room  temperature  and  monitored  the 
absorption  edge.  Their  result  was  the  relation 

Eg  »  1.45  ♦  0.0126P  -  3.77’‘10‘5p2  (eV).  (70) 

where  P  is  in  kbar. 

Another  external  variable  affecting  Eg  is  temperature.  Varshni^23) 
proposed  a  temperature  dependence  given  by 

Epd) »  Ei.(0)  -  ocT2/(T*^)  (eV).  (70 

with  «  =  8.87I>‘10''^  eV/K,  and  ^  =  572  K. 

Additional  results  were  then  considered  and  the  values  of  the  constants 
were  suggested  to  be  «  =  5.405*10”^  eV/K  and  ^  =  204  K. 


V.  nAGNETIC  FIELD  EFFECT 


A.  Simple  Bands 

As  an  introduction  to  the  effect  of  a  magnetic  field  on  the  absorption 
coefficient  we  first  consider  a  simple  case:  a  free  particle  of  mass  m  and 
charge  -e  and  spherical  non-degenerate  bands,  with  the  magnetic  field 
uniform  and  in  the  *z  directioa^^^  Using  the  Landau  gauge,  the  vector 
potential  is  fl  =  H(0x0),  and  the  Hamiltonian  is 

K  =  (2m)"‘lp-(e/c)RP 
=  (2m)‘‘[p2^*p^2*(Py-eHx/c)2L 

Using  the  definition  of  the  cyclotron  frequency,  w^-eB/mc,  the 
Hamiltonian  becomes 

H  =  (2mr 

and  the  Schroedinger  Equation,  H'f  =  Ef,  can  be  written  as 

or. 

(tsV2mX92t(x)/9x2>»(  li2j(2^2m)t(xM2m)  "Kfiky^mo>j^)2t(x)  =  E’Kx) 


where  Py,2«*ky  2- 


Rearranging  this  expression. 


-(^\V2mX92♦(x)/^x2)■»[(fn6>  ^/2Xx**ICg/m<0c)^ 

Hi2k22/2mlt(x)  s  Et(x).  (72) 

Equation  (72)  has  the  form  of  a  simple  harmonic  oscillator  equation  in 
one  dimension,  with  equilibrium  position  at  x‘  -  fiky/mo>c  >  -cfiky/eB.  and 
natural  frequency  6)^  The  eigenvalues  are 

E„  =  (nH/2)fi<i)c*ftk22/2m. 

Since  the  equilibrium  position  of  the  oscillator  must  be  within  the 
material  (consider  a  parallelepiped  of  dimensions  L^.  Ly,  L2).  a  limitation 

is  imposed  on  the  range  of  ky,  i.e.  -1^/2  <  x'  <  L^/2,  or  upon  writing  x*  in 


terms  of  ky, 


-(Lx/2XeH/oh)  <  ky  <  (Lx/2)(eB/cri). 


There  are  l-y/2TT  allowed  values  of  ky  and  we  have  a  degeneracy  of 


levels,  for  fixed  k2,  of  (Lj,Ly/27TXeB/ch).  If  spin  is  neglected,  this 
expression  must  be  multiplied  by  a  factor  2.  Further,  the  number  of 


allowed  values  of  k^.  such  that  |k2|<  k20.  is  L^k^o/Tf  which  also  gives  the 

allowed  values  of  Ej^  <  Ej^o  =  h^k2q^/2m.  Solving  for  k^o  and  substituting, 
we  obtain  for  the  number  of  allowed  values 

Writing  E|(  for  Ej^o  and  differentiating  with  respect  to  results  in  the 
number  of  allowed  energy  levels  in  dE|^.  Thus. 

N(Ek)dE,j  =  (L2/27m)(2m)*/2(E,j)‘»/2  (73) 

in  order  to  find  the  distribution  of  ail  the  levels,  we  sum  over  the 
distribution  functions  and  use  the  degeneracy  factor.  The  density  of  states 
then  becomes,  for  simpie  bands,  using  An  =  0. 

Jcv=(L*Lj|/2irXeH/oliXLz(2mr)'/V2milZ,^l.<i)-Eg-(n‘l/2)  ftWo-r A  (74) 

•s 

Where  m^  is  the  reduced  mass,  and  *  a>y  ^eB/cny  . 

Equation  (74)  is  generally  written  as 

Jcv  =  (V/47T2)(eH/cft)l(2mrV''2/ftl2nlftw-Eg-(n+l/2)  (75) 

since  volume. 

Using  Eq.  (75),  the  absorption  coefficient  becomes 


o<(<i>,H)=(47rV/ncm2<DV)  \  K'  |  ^*(V/4Tr2XeB/chX(2mr)’^Vftl 


or,  upon  simplifying. 


o<(<l>,H)=le2Ep/6ncm(^S«)KcH/cl^X2my^)’^2^^f <i)-Eg-(n*l/2)  (75) 

Using  the  definition  of  (0^^  *  eH/cny.  we  write 

oc(<o.H)  =  le2Ep(2mr)3/V6ncm(ft<i>)ft2Xiio>y2)ErJ[f  (i)-Eg 

-(n»i/2)^i<i)crrA  (77) 

Equation  (77)  compares  with  the  zero-field,  no  spin  case  given  by 

o<o(a>)=lc^(2mr)^V6ncm(fi<i>)ft2Kti<«>-Eg)’^2, 
or,  upon  writing  o<(<i),0)  =  C(ft<i)-Eg)'^2,  we  see  that 

«(<i).H)  =  (C/2)li<i>cr<fla>-En)‘'/2, 

where  Ep  «  Eg  ♦(n+l/2)fi<i)^;r  • 

A  dimensional  check  shows  that 

o<((i>.B)  [e2Epmr^/Vcm(ft<o)h2Kft(Dcr/E>/2) 

-V  [ExL»cExE3/2x(t3/L3)xE  /  (L/t)x(E/L2)xt2xExE2xt2xE»/2l  -v  L"’, 


Examination  of  Eq.  (77)  shows  that  the  application  of  a  magnetic  field 
changes  the  energy  bands  from  the  simple  situation  shown  in  Fig.  la,  to  the 
more  complicated  case  shown  in  Fig.  1b.  Further,  using  Eq.  (77).  we  can 

find  o<((i>,B).  Using  the  values  for  m^  and  m^  determined  earlier 
(m^.«0.067m.  myS!0.505m,  fTy«0.059m)  and,  for  purposes  of  illustration 
choosing  a  field  of  l  Weber/m^  =  lO^  6.  we  have,  using  *6)^  =  ft(eB/cm). 

1.158*IO’'<(B/(mr/m)l  (eV)  =  1.963>‘10"3  (eV). 

Ihea  for  example,  using  the  value  of  just  calculated,  we  have 
o<(«.l)M23/2c2Epm’''20.0593''2/6noR2Xf<«>cr/2X<i(i))‘’ 

*ZrjHi(i)-Eg-(n*1/2)fia>^.rr’''2 

«(Z163*10«X0.014X9.82*10''’)(<ia>)*’B<^(i)-Eg-(ml/2)<^<OcrP/2 

or,  o<(<a.l)«29.635(f  <i>)'’2n[ft<o-Eg-(n^I/2)  ftCDcrr’^^.  (78) 


The  energy  levels  of  degenerate  bands  in  the  absence  of  external 
perturbation  were  calculated  using  Kane’s  k  .  p  method.  Luttinger  and 
Kohrf2^)  Investigated  the  problem  of  cyclotron  resonance  for  degenerate 


bands,  and  Luttinger  described  a  general  method  for  determining  the 
energy  levels  of  holes  in  a  degenerate  band.  Since  in  our  samples  the 
degeneracy  was  lifted  by  strain  (see  Sect.  Vl-F),  we  were  essentially 
dealing  with  the  simple  band  situation  described  above. 

B.  Eicitons  in  Hagnetic  Fields 

Group  ni-V  compounds  generally  have  large  dielectric  constants  and 
comparatively  small  effective  masses  at  the  conduction  and  and  valence 
band  extrema.  Therefore  excitons  in  these  materials  have  a  large  spatial 
extent,  with  the  exciton  radius  given  by 

where  p  Is  the  reduced  effective  mass,  ^  is  the  dielectric  constant.  The 
ground  state  energy  is  below  the  gap  energy  by  the  amount. 

=e'^p/2«o2fi2 

As  noted  earlier,  the  exciton  effect  is  mainly  apparent  in  a 
modification  of  the  absorption  edge,  and  spectra  do  not  generally  show 
much  structure.  Thus,  although  the  valence  band  maximum  In  GaAs  is 
complex,  the  absence  of  detail  in  the  spectra  means  that  a  simple 
hydrogenic  model  is  sufficient  to  obtain  the  exciton  states  and 


energies^^sj  usually  spectra  show  only  transitions  associated  with  the 
lowest  exciton  state,  and  interband  magneto-optical  phenomena  are 
basically  explained  in  terms  of  the  Landau  levels,  possibly  perturbed  by  the 
electron-hole  Interaction/^^^ 

In  quantum  well  structures  the  exciton  can  be  considered  to  be 
two-dimensional  with  a  binding  energy  approximately  four  times  larger 
than  in  the  three-dimensional  case.^2?)  Hamiltonian  can  be  written 

35(28) 

H  =  Ke  ♦  Kh  ♦  Kx. 

where  Kg  and  Hf,  represent  electron  and  hole  motion  in  the  2-dircction,  and 

Hx  =  PV2M  ♦  pV2|i  -e2/«(r2*z2)’''2  (79) 

In  Eq.  (79).  P  =  pg+p^.  M  =  mg+m^.  ji  is  the  reduced  mass,  and  r  =  rg- 

2  =  Zg-  Zf,.  Naturally  M  and  k  have  different  values  depending  on  whether 
heavy  or  light  holes  are  being  considered,  p  is  the  momentum  conjugate  to 


r,  and  generally  the  total  momentum.  P,  is  taken  to  equal  zero.  In  order  to 
include  the  external  magnetic  field,  we  can  rewrite  Eq.  (79)  as 


I 


H  =(2mgr’[pe-(e/2c5  ♦  (2m^P[Ph*(e/2c)B 

-e2/e|rg-rj^|.  (80) 

When  the  magnetic  field  is  strong  enough  to  allow  omitting  the 
Coulomb  term  in  Eq.  (80),  we  are  left  with  an  expression  leading  to  the 
Landau  levels.  This  was  considered  in  the  previous  sectioa  The 
two-dimensional  situation  (well  width  negligible),  for  unrestricted  values 

of  the  magnetic  field  (B^).  was  investigated  by  Akimoto  and  Hasegawa.^^^ 

They  determined  approximate  expressions  for  the  energy  as  a  function  of 
in  the  ’low'  and  "high'  field  regimes.  Low  field  has  been  loosely  defined 


as  the  region  where 


<i(eB/2MC)  «  Me^/2f»*eto2 


and  high  field  where  the  converse  holds.^26)  for  -|ow"  field 

Erj/Ry*'=-(ml/2)"24(5/8)(ml/2)‘^(li<«>o)2/(Ry‘')2  ♦...  .  n=0.l.2 .  (81) 

and,  for  the  "high"  field. 

E^y*‘=  (2ml)t»<i>o/Ry*‘-3lf(Oo/(2n-l)Ry‘'r/2*.. ,  n=0.1.2 (82) 

In  these  two  expressions  the  following  definitions  apply: 


*0)0  **(eB/2pc).  and  Ry“  *  |ieV2fi2cp2 

For  GaAs  we  remain  in  the  low  field  regime  for  H<60kG,  and  since  we  are 
concerned  with  the  lowest  level  exciton.  we  can  rewrite  Eq.  (81)  as 

Eo  «  MRy**  ^  (5/8Xl/2)^(*o>o)VRy^  (83) 

in  Fig.  2  we  illustrate  the  magnetic  field  dependence  of  energy  in  the 
case  of  Landau  levels  (simple  bands)  using  the  relation  E^  =  (n»l/2)*eB/]ic. 


and  in  the  case  of  excitons  using  Eq.  (83). 


1/2)liw 


dependence  for  Landau  energies  (top)  and  excitons 
lerqu  scales  are  different  to  show  effect. 


VI.  EXPERinENT 


A.  Experimental  Objective 

The  purpose  of  this  experiment  was  to  develop  a  method  to  allow  the 
optical  investigation,  at  low  temperatures,  of  semiconductor  structures 
subjected  simultaneously  to  high  pressure  and  high  magnetic  field.  This 
objective  was  met  by  the  experiment  described  in  the  following  paragraphs. 
Although  we  demonstrated  the  feasibility  of  the  experimental  method  by 
performing  laser  induced  luminescence  measurements,  the  experiment  is 
readily  adaptable  to  investigations  involving  absorption,  transmission,  and 
excitation  spectroscopy. 

B.  General  Experiment  Description 

The  novel  feature  of  this  experiment  is  that  optical  access  to  the 
sample  to  be  investigated  is  by  means  of  an  optical  fiber.  This  alleviates 
many  of  the  problems  often  associated  with  measurements  involving  high 
pressures.  e.g.  limitations  imposed  by  pressure  chamber  and/or  Dewar 
windows,  transmissivity  of  the  pressure  medium,  and  also,  in  the 
configuration  we  used,  it  allows  simultaneous  determination  of  magnetic 
and  pressure  effects. 

In  our  experiment  the  sample  was  attached  directly  to  one  end  of  the 


optical  fiber  and  suspended  in  the  pressure  chamber.  Fiber  access  to  the 
pressure  domain  was  through  a  pressure  plug.  The  chamber,  at  the  end  of  a 
probe,  was  inserted  into  the  magnetic  field  on  the  axis  of  a 
superconducting  magnet  (solenoid).  The  sample  was  illuminated  by  an 
Argon^  Ion  laser  beam  which  was  coupled  into  the  other  end  of  the  fiber 
through  a  microscope  objective.  The  luminescence  signal  from  the  sample 
returned  through  the  same  fiber,  and  the  collimated  signal  was  split  off  and 
directed  to  a  monochromator  and  detector.^^o) 

The  magnet  is  capable  of  providing  fields  to  65  kG.  and  the  pressure 
system  was  used  to  4  kbar  (*^60000  psi).  Hydrostatic  pressure  is  exerted 
on  the  sample  by  growing  a  uniform  crystal  of  helium  around  it  under 
controlled  conditions. 

C.  Preliminary  Investigations 

Certain  preliminary  work  was  required  to  test  the  feasibility  of  the 
proposed  method.  This  included  determination  of  fiber  suitability  at 
temperatures  <  4  K.  and  preparation  of  fiber  ends  (polishing).  Additionally 
a  means  of  attaching  the  sample  to  the  fiber  had  to  be  found,  as  well  as 
providing  a  seal  around  the  fiber  where  it  entered  the  pressure  regioa  That 
optical  losses  had  to  be  kept  to  a  minimum  was  understood. 


Fiber  eixls  were  polished  to  0.3  pm  smoothness  using  progressively 
finer  grit  abrasives,  with  the  acceptance  test  being  visual  smoothness 
under  lOOx  magnification.  A  Gaussian  beam  propagation  calculation  was 
used  to  select  the  objective  to  match  the  beam  waist  to  the  fiber  core 
diameter.  Positional  adjustments  to  the  input  end  of  the  fiber  were  made 
by  means  of  a  3*degree  freedom  micrometer  stage  holding  the  fiber,  and 
transmission  was  maximized  initially  by  using  a  power  meter  at  the  sample 
end  of  the  fiber.  A  beam  splitter  was  used  along  with  folding  mirrors  (Ftl) 
to  direct  the  return  signal  from  the  sample  (initially  a  mirror)  to  a 
spectrometer.  This  Inital  arrangement  is  shown  in  Fig.  3. 

In  this  arrangement.  Fhs  1*3  were  used  to  position  the  beam  on  the 
microscope  objective  which,  in  tura  placed  the  beam  waist  on  the  the  end 
of  the  fiber.  The  return  signal  placement  was  adjusted  by  means  of  mirrors 
Fn4  and  Fn5.  The  lens.  L.  focused  the  signal  on  the  spectrometer  input  slit. 
After  initial  positioning  of  the  optical  elements  resulted  in  a  visible  spot 
on  the  screen,  assuming  the  spectrometer  was  properly  adjusted,  the  signal 
was  maximized  (visually)  by  adjustment  of  fiber  position  and  the  beam 
splitter.  BS.  Greater  alignment  accuracy  was  then  obtained  by  replacing  the 
screen  with  a  photomultiplier  tube  (PMT).  Displaying  the  PUT  output 


on  an  oscilloscope,  fine  adjustments  were  compieted.  For  this  step  a 
chopper  wheel  was  placed  between  FI12  and  Fn3.  and  the  beam  stop  was 
replaced  by  a  photodiode  to  provide  the  external  trigger  signal  to  the 
scope.  This  set-up  performed  satisfactorily  with  a  CdS  sample 

The  laser  in  this  preliminary  test  was  a  Spectra  Physics  Model  164 
Argon^lon  laser  with  Model  265  Exciter,  used  to  pump  a  Spectra  Physics 
Model  375  dye  laser  with  Rhodamine  66  dye.  The  spectrometer  was  a 
Jobin-Yvon  Model  H-20V.  and  the  PMT  was  an  RCA  7767  with  S-1  response. 

Although  the  set-up  shown  in  Fig.  3  worked  well,  two  significant 
problems  soon  became  apparent.  Sample  adhesion  to  the  fiber  was  not 
good,  and  signal  losses  when  the  end  of  the  fiber  and  sample  were 
immersed  in  liquid  nitrogen  were  noticeable.  The  sample  adhesion  problem 
was  solved,  after  experimentation,  by  thinning  GE  Adhesive  and  Insulating 
Varnish  No.  7031  with  a  50/50  mixture  of  toluene  and  methanol.  A  general 
guide  required  I  part  varnish  and  5-10  parts  thinner  by  volume.  Optimum 
results  were  obtained  when  the  mixture  reached  a  consistency  to  form  a 
bead  at  the  end  of  a  pin  without  dropping  off.  The  adhesive  then  had  to  be 
used  to  be  used  promptly  since  it  dried  rather  quickly.  The  magnitude 
of  the  signal  loss  problem  did  not  become  apparent  until,  in  the  final 


experimental  configuration,  approximately  four  feet  of  the  fiber  (inside  the 
probe)  were  cooled  to  liquid  nitrogen  temperature  and  below. 

This  particular  fiber  (SS-600.  General  Fiber  Optics.  Inc.)  supposedly 
consisted  of  a  fused  silica  core  wih)  doped  silica  cladding,  and  the  indices 
were  supposed  to  trade  each  other  with  temperature.  We  therefore  felt 
initially  that  the  poor  or  nonexistent  fluorescence  signals  were 
attributable  to  other  causes.  e.g.  sample  quality  and/or  adhesion  and 
placement  of  the  sample  on  the  fiber.  A  test,  which  should  have  been 
performed  earlier,  showed  that  transmitted  power  became  negligible  when 
approximately  four  feet  of  colled  fiber  where  immersed  in  liquid  nitrogen 
Our  conclusion  is  that,  in  place  of  the  doped  silica  cladding  advertised,  our 
fiber  had  a  polymer  cladding  whose  index  had  changed  sufficiently  at  76  K 
to  allow  most  of  the  light  to  escape  through  the  sides.  An  Aheoated 
quartz  fiber  manufactured  by  Hughes  Aircraft  Company,  designated  P-27, 
with  a  core  diameter  of  100  pm  solved  the  problem.  This  fiber  showed 
negligible  loss  at  76  K  and  was  used  henceforth.  A  further  advantage  of 
P-27  was  that,  upon  cleaving,  it  presented  a  smooth  surface  without 
polishing.  The  Al  coating  was  cleared  from  the  end  surface  using  a 
laboratory-grade  NaOH  solution 


D.  Experimental  Configuration  and  Equipment 

Since  the  optical  system  used  in  the  preliminary  tests  performed  well, 
it  was  used  with  only  minor  ad^tatic^  for  the  experiment.  The  changes 
were  required  because  the  positions  of  the  laser  and  monochromator  used 
for  the  data  runs  were  dictated  by  available  bench  space.  Figure  4  shows  a 
diagram  of  the  optics  configuration  In  this  experimental  configuration,  all 
optical  elements  beyond  FI13  were  motmted  on  a  platform  adjusted  to  the 
height  of  the  monochromator  input  slit.  Thus  FMs  1  -  3  are  used  to  raise 
the  laser  beam  to  the  appropriate  level,  and  only  lateral  adjustments  to  the 
other  optical  elements  are  required.  During  all  data  runs  the  Argon^  Ion 
laser  was  tuned  to  5145  angstrom  (2409.9  meV)  and  an  output  power  of  1.2 
-  1.4  w  (cw).  The  beam  splitter  reflects  4X  of  the  incident  light.  The  red 
filter  with  a  cu^  ^ff  near  6000  angstrom  prevented  reflected  laser  light 
from  entering  the  SPEX  0.75  m  double  monochromator.  A  RCA  7102 
phot')detector  was  used  at  the  output  slit  with  a  Fluke  415B  HV  Power 
Supply,  and  its  output  was  recorded  on  a  HP  7044B  x-y  Recorder. 

The  fiber  was  then  sealed  in  the  pressure  plug  and  inserted  into  the  probe. 
Both  components  were  made  of  heat-treated  BeCu.  A  sketch  of  the  plug  is 


Fig.  4.  Placement  of  optical  elements.  Laser  beam  Is  coupled  Into  the 
optical  fiber  to  Illuminate  the  sample,  and  the  return  signal  Is  directed  to 
the  monochromator. 


shown  in  Fig.  5,  with  relevant  dimensions  givea  The  probe  has  an  outside 
diameter  of  0.250  inch,  and  and  inside  diameter  of  0.0625  inch.  The 
diameter  of  the  chamber  opening  is  0.125  inch. 

The  seal  between  the  probe  shaft  and  the  bomb  Is  ensured  by  exact 
machining  of  angles  and  compression  provided  by  threaded  fittings.  The 
fiber  was  sealed  in  the  pressure  plug  with  epoxy.  The  magnet  is  contained 
in  a  triple  Dewar,  as  shown  in  Fig.  6.  and  driven  by  a  Varlan  V4102 
Superconducting  Magnet  Power  Supply.  Level  sensors  in  the  outer  and 
middle  Dewars  and  a  pressure  sensor  in  the  inner  Dewar  provided  necessary 
informatioa 

Hydrostatic  pressure  on  the  sample  was  provided  by  a  method  dubbed 
'Solid  Helium  Pressure  Generation  Technique'.^^’^  Basically  this  technique 
involves  maintaining  pressure  on  liquid  helium  as  it  solidifies  slowly  from 
the  bottom  of  the  pressure  chamber  and  surrounds  the  sample.  After  the 
solid-liquid  interface  has  passed  the  sample,  the  probe  can  be  lowered  to 
the  appropriate  position  since  additional  cooling  can  be  considered  to 
occur  at  constant  volume  due  to  blockage  of  the  probe  capillary.  The 
pressure  system  Is  shown  in  Fig.  7.  The  system  was  charged  initially  to 
'^2500  psi  with  helium  gas.  Pump  I  was  then  used  to  a  pressure 


Fig.  5.  The  pressure  plug  (BeCu),  encasing  the  optical  fiber,  provides  the 
seal  between  the  high-pressure  region  and  the  ambient  atmosphere  and 
allows  optical  access  to  the  pressure  regioa 
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pressure  generating  system. 


of  '^20,000  psi,  and  thereafter  maximum  pressure  was  attained  using  Pump 
2.  Detailed  operating  instructions  are  included  under  Experimental 
Procedures.  After  the  system  had  been  prepared  for  a  data  run,  the  actual 
measurement  and  data  recording  was  controlled  by  a  Computer  Automated 
Measurement  and  Control  (CAMAC)  system  developed  by  Dr.  E.  D.  Jones  of 
Sandia  National  Laboratory.^^^^ 

E.  Experimental  Procedures 

I.  Fiber  Preparation 

A  length  of  fiber  was  cut  from  the  roll.  The  length  had  to  be 
sufficient  to  allow  raising  the  probe  vertically  above  the  Dewar  assembly 
while  the  other  end  of  the  fiber  was  mounted  in  the  adjustable  fiber  holder. 
Both  ends  of  the  fiber  were  cut  using  a  fiber  scribe,  and  examined  under  a 
IOOj*  microscope  for  flatness.  If  carefully  done,  fiber  P-27  did  not  need  to 
be  polished.  The  jagged  aluminum  coating  was  then  removed  from  the  ends 
with  a  NaOH  solutioa  The  fiber  was  inserted  through  the  pressure  plug  and 
sealed.  The  seal  was  provided  by  Eccobond  285  Black  with  Catalyst  24LV. 
mixed  in  a  ratio,  by  weight,  of  100  parts  Eccobond  to  7-8  parts  of  Catalyst 
24LV.  Curing  time  at  room  temperature  Is  ~24  hours.  In  order  to  locate 
the  sealant  in  the  critical  region  in  the  narrow  part  of  the  pressure 


plug,  a  configuration  as  shown  in  Fig.  6  was  used.  The  epoxy  was  inserted 
into  the  sleeve  attached  to  the  nipple  of  the  pressure  plug  by  using  a 
pipette  and  was  pulled  down  into  the  plug  by  a  vacuum.  The  back  seal  during 
this  process  was  provided  by  a  rubber  cap  through  which  the  fiber  had  been 
fed  with  a  syringe  which  was  then  withdrawn.  The  vacuum  was  left  on  the 
system  for  approximately  ten  minutes  after  which  the  epoxy-filled  upper 
sleeve  was  carefully  cut  off  and  the  fiber  cleaned.  This  had  to  be  done 
since  the  plug  does  ftot  fit  into  the  pressure  joint  with  the  sleeve  attached. 
The  assembly  was  then  left  to  cure. 

2.  Sample  Attachment 

After  curing,  the  sleeve  was  removed  and  the  fiber  with  pressure 
plug  was  inserted  into  the  probe  and  secured  in  the  pressure  joint  with  the 
proper  fitting.  The  compression  exerted  when  the  fitting  was  tightened 
served  to  improve  the  seal  around  the  fiber.  If  the  positon  of  the  plug  on 
the  fiber  was  carefully  measured,  slightly  over  one  inch  of  fiber  protruded 
from  the  probe.  The  probe  was  then  mounted  vertically  above  an  elevator 
stage.  The  end  of  the  fiber  and  the  sample  were  cleaned  with  methanol  (and 
lens  paper)  and  the  sample,  positioned  on  the  stage,  was  brought  up  into 
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Fig.  8.  Pressure  plug  holder  used  during  sealing  around  the  optical  fiber 


contact  with  the  fiber.  Since  the  sample  size  was  usually  about  1/16 
inch  square,  we  generally  tried  to  center  the  fiber  on  the  sample.  When  the 
two  were  in  contact  the  adhesive  was  prepared  (GE  7031  plus  501K 
toluene/50%  methanol).  Upon  reaching  the  desired  consistency,  it  was 
applied  to  two  sides  of  the  fiber  at  the  juncture  with  the  sample.  Care  had 
to  be  taken  to  ensure  that  the  sample  did  not  become  attached  to  the 
platform,  or  that  the  fiber  shifted  on  the  sample.  If  this  occurs,  the 
adhesive  must  be  removed  (acetone)  and  fiber  and  sample  must  be  recleaned 
before  starting  anew.  The  adhesive,  when  properly  prepared,  dries  quickly 
and  allows  attachment  of  the  pressure  chamber  within  one  hour, 

3.  Alignment 

The  probe  was  then  ready  for  insertion  into  the  Dewar.  After  this 
had  been  accomplished,  the  other  end  of  the  fiber  was  positioned  in  the 
fiber  holder.  The  Dewars  could  be  filled  at  this  point  if  only  magnetic  field 
measurements  were  planned,  but  a  different  procedure  applies  for  pressure 
measurements.  With  the  laser  operating,  the  monochromator  was  scanned 
to  locate  the  luminescence  peak,  and  adjustments  of  the  fiber  position  and 
optical  elements  were  used  to  maximize  the  signal.  If  necessary, 
attenuators  could  be  placed  in  the  beam  path. 
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4.  nsasuTBiiwnts  of  riagnetic  Field  Effects 


With  the  system  prepared  as  outlined  above,  the  magnet  was 


adjusted  to  the  desired  field  strength  and  the  energy  interval  to  be  recorded 


was  selected.  The  monochromator  was  set  to  the  starting  point  of  the 


scan  interval  and  data  recording  was  controlled  by  the  CAMAC  ^tem. 


S.  tleasiirements  of  Pressure  Effects 


is 
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When  pressure  measurements  were  intended,  the  procedure  differed 
in  that  probe  insertion  had  to  include  specific  steps.  To  start  with,  the 
Dewars  had  to  be  filled  to  the  appropriate  levels  before  insertion  was 
begua  Additionally,  the  system  must  be  pressurized  and  a  Cu-constantan 
thermocouple  had  to  be  attached  to  the  chamber  to  give  an  indication  of 
temperature.  Step-by-step  Instructions  for  achieving  the  desired  pressure 
with  the  system  shown  in  Fig.  7  are  given  below.  A  full  bottle  of  helium 
gas  (^^2500  psi)  was  connected  to  the  system.  Referring  to  Fig  7, 

a.  Open  V2  and  Vt  Bleed  (1/2  turn) 

b.  Open  VI  Charge  until  hissing  sound 

c.  Adjust  hiss  for  slow  rate  with  VI  Charge 

d.  Close  VI  Bleed 

e.  Wait  for  system  to  charge  to  bottle  pressure 


f.  Close  V2  and  VI  Charge 

g.  Open  VI  Bleed 

h.  Close  valve  on  Pump  t 
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i.  OpcnV3 

j.  Use  Pmp  1  to  '^20,000psi  (or  operating  pressure,  if  less) 

k.  Close  V3  and  open  valve  on  Pump  i 

l.  Open  V4  Charge 

m  Use  Pump  2  to  desire  maximum  pressure. 

The  system  was  depressurized  by  following  the  steps  listed  below, 
after  closing  the  valve  on  the  helium  bottle. 

a.  open  V4  Bleed  SLOWLY  (to  >^20,000  psi) 

b.  Close  V4  Bleed 


c.  CAREFULLY  open  V3  for  slow  decrease  to  ~2500  psi 


d.  Close  V3 


e.  Close  VI  Bleed 


f.  OpenV2 

g.  Open  VI  Bleed  until  slow  hiss. 

With  the  system  pressurized  according  to  the  outlined  procedure,  the 
probe  was  lowered  slowly  into  the  Dewar,  while  monitoring  the 


temperature  with  the  thermocouple.  It  is  important  to  lower  the  probe 
slowly  enough  that  the  thermocouple  and  pressure  chamber  can  be 
considered  to  be  in  thermal  equilibrium.  When  the  freezing  temperature  of 
helium  at  the  selected  pressure  was  approached,  the  rate  of  descent  was 
slowed  to  approximately  1  inch/hour  to  allow  a  uniform  helium  crystal  to 
grow  around  the  sample.  Once  the  solid-liquid  interface  passed  the  sample, 
the  rate  of  descent  was  increased.  During  this  entire  process  the  probe 
was  suspended  from  a  pulley  arrangement,  and  a  plastic  bag  was  used 
around  the  probe  at  the  top  of  the  Dewar  to  prevent  the  access  of  air.  The 
thermocouple  was  referenced  to  liquid  nitrogen  and  the  freezing 
temperature  of  helium  was  taken  from  P  vs  T  charts. 

Once  the  probe  was  in  position,  alignment  and  recording  procedures 
were  as  before.  In  order  to  perform  a  measurement  at  another  pressure, 
the  probe  was  raised  until  the  temperature  of  the  chamber  was  above  the 
helium  melting  point,  the  pressure  was  adjusted,  and  the  process  was 
repeated.  This  method  is  described  in  detail  in  Reference  (31). 

F.  Samples 

A  large  number  of  samples  were  examined  in  literally  hundreds  of 
experimental  runs  during  the  course  of  this  experiment.  Naturally  many  of 
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these  runs  were  required  to  establish  experimental  procedures  during  the 
initial  stages  of  the  experiment  and  can  not  be  counted  as  actual  data  runs. 
The  samples  ranged  from  bulk  materials  (GaAs,  InGaAs)  to  single  and 
multiple  quantum  well  structures,  as  well  as  strained-layer  superlattices. 
The  quantum  well  and  super  lattice  structures  were  constructed  of 
alternating  layers  of  different  semiconductor  materials  with  various  layer 
widths  and  numbers  of  layers.  Samples  were  predominantly  of  the  type 

GaAs/AIAsorGaAs/ALj,Gai.jjAs.  and  lr^Ga,-^As/GaAs.  where  the  subscript 

X  refers  to  the  percentage  of  the  particular  component  present  in  the  layer. 
The  number  of  layers  ranged  from  the  case  of  a  single  layer  of  GaAs 
sandwiched  between  two  layers  of  AlAs  in  the  single  quantum  well  to 
multiple  wells  of  six  to  twelve  periods,  and  to  superiattices  with  fifty  or 
more  periods.  The  typical  sample  included  a  substratum,  a  buffer  layer,  the 
alternating  well  and  barrier  layers,  and  a  cap  layer. 

L,  R.  Dawson  (SNL)  grew  the  structures  by  molecular  beam  epitaxy  on  a 
<100>  GaAs  substratum.^55'5^^  The  substratum  was  placed  In  a  pressure- 
and  temperature-controlled  growth  chamber,  and  effusion  cell  shutters 
were  opened  for  specific  time  intervals  to  control  the  molecular  beam 


depositions  for  the  different  layers.  The  advantage  of  this  computer 
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controlled  process  lies  in  the  accurate  control  it  affords  over  layer 
thickness  and  composition.^^®*^®^ 

The  potential  energy  profile  of  such  a  heterostructure  is  shown  in  Fig. 

9.  In  the  case  of  6aAs/AI^6ai.^As,  the  6aAs  wells  are  bounded  by  AIGaAs 

barriers  whose  height  is  determined  by  the  fraction  of  Al  present.  For 

values  of  X  <  0.4,  Alx6a|.xAs  maintains  a  direct  gap  at  Electron  and 

hole  energies  are  quantized,  with  levels  for  "heavy*  and  "light*  holes  in  the 
valence  band  wells,  and  transitions  between  electron  and  hole  states 
conform  to  the  selection  rule  An  =  0.  In  this  particular  structure  strain 
effects  are  negligible  since  a  good  lattice  match  exists.  The  wells  can  be 

quite  rectangular  with  nominal  depth  of  of  AEg  for  the  conduction 

band,  and  of  AEg  for  the  valence  band,^’®^  and  the  gap  at  T  =  0  K  given 
by  (39) 

Eg  *  t.52  ♦  t.52X  (eV).  (84) 


The  band  gap  difference.  AEg.  has  been  determined  empirically  to  be^^r) 


AEg  =  1.155X  ♦  0.37x2  (eV). 


(85) 


In  the  case  of  ln^Gai_^As/GaAs  structures,  the  GaAs  layers  form  the 
barriers  for  the  InGaAs  wells,  where  the  amount  of  In  present  determines 


the  minimum  gap  energy.  In  this  situation  GaAs  remains  transparent  near 
the  InGaAs  gap.  and  this  gap  is  given  by 

Eg  «  1.52  -  1.601X  ♦  0.536x2  (eV).  (86) 

These  structures,  unlike  GaAs/AIGaAs.  are  not  lattice  matched. 
Although  this  can  result  in  dislocations  in  thick  layers,  the  mismatch  is 
fully  accomodated  by  strain  if  the  layers  are  thin  enough.  If  the 
orientation  of  the  superlattice  is  in  the  <100>  direction,  the  strain  can  be 
decomposed  into  hydrostatic  and  uniaxial  <I00>  components.^^®^  The 

hydrostatic  strain  results  in  an  increase  in  Eg.  while  the  uniaxial 

component  lifts  the  degeneracy  of  the  valence  band  maxima.^^’^  The 

expression  for  Eg.  Eg.  (86).  does  not  include  this  effect.  A  value  of  x  =  0.25 

corresponds  to  a  lattice  mismatch  of  1.8*.^^°^  and  the  value  of  Eg  as  a 

function  of  lattice  constant  is  shown  in  Fig.  10. 

In  simple  photoiuminescence  measurements,  these  samples  are 
expected  to  show  a  single  peak  from  the  InGaAs  well  due  to  diffusion  of 
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Fig.  10.  Band  gap  of  GaAs-lnj^6a|_j^  As  as  a  function  of  lattice  constant 

parallel  to  the  interface.  The  curve  labeled  Alloy  refers  to  bulk  InGaAs,  and 
the  dimensions  listed  apply  to  the  layers  of  both  materials  (in  different 
samples). 
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the  carriers  to  the  well  before  recombinatioa  The  transition  involves  the 
first  quantum  well  levels  and  may  also  Involve  exciton  effects.^33.42) 
position  of  the  peak  energy  is  a  function  of  the  amount  of  In  and  the  well 
width.  The  application  of  hydrostatic  pressure,  as  in  our  experiment,  is 
then  expected  to  shift  this  peak  to  higher  energy  without  altering  its  shape. 
For  GaAs/AIGaAs  structures,  experiments  have  shown  that  in  undoped 
samples  exciton  recombination  is  the  dominant  radiative  process,  while 
intrinsic  band-to-band  recombination  dominates  in  modulation  doped 
samples. 

A  major  difference  between  the  multilayer  heterostructures  and  bulk 
alloy  materials  is  the  effect  on  optical  properties  resulting  from  the 
periodicity  of  the  layers,  especially  when  the  layer  thickness  approaches 
the  DeBroglie  wavelength  of  the  electronic  particle.  This  limitation  on 
carrier  motion  in  the  direction  normal  to  the  layers  is  then  analogous  to 
the  confinement  of  carriers  to  a  one-dimensional  potential  well.  If  the 
well  depth  is  considered  to  be  infinite,  then  the  energy  can  be  written,  for 

Ljj,  Lg,  as 

E  =  Er,  ♦  (^i2/2m)(kx2^ky2).  (87) 


where  n  =  1. 2, 3 . 

For  each  value  of  a  two-dimensional  energy  band  forms  in  the  k^ky  plane 

and.  as  decreases,  the  energy  of  the  discrete  band  states  will  increase. 

This  is  the  quantum  size  effect.  The  confinement  effect  on  excitons  has 
been  discussed  earlier. 

When  the  barrier  thickness  is  large,  optical  effects,  e.g.  absorption 
are  basically  described  by  the  single  well.  However,  when  wells  are 
brought  close  together  by  reducing  the  thickness  of  the  finite  barriers,  the 
coupled- well  regime  is  reached.  The  interaction  between  the  bound  states 
then  splits  these  levels.  For  a  small  number  of  coupled  wells  the  splitting 
can  be  detected.  For  example,  transitions  involving  symmetric  and 
anti-symmetric  combinations  of  the  n  =  1  bound  states  yield  a  fine 
structure  within  this  level.  As  the  number  of  coupled  wells  increases,  this 
discrete  structure  tends  to  disappear  and  a  broad  band  (several  meV) 
appears  in  the  superlattice  regime. 

We  now  calculate  the  band  edge  effective  masses  for  ln^Gat_^As.  In 

order  to  do  this  we  used  the  expressions  for  energies  developed  by  the  k'p 
method,  as  given  in  Eqs.  (37).  and  first  obtain  the  vaiues  of  the  constants 


for  this  materials.  Interpolating  between  the  E(r|)  values  of  InAs  and 
GaAs/^^^  we  have,  for  x  =  0.2: 

E(r,5v)  =  0.  E(r,c)  =  1.28,  E(ri5c)  =  =  ’2.  fp  =  25.2.  a  =  0.35. 

all  in  eV.  Using  these  values,  we  obtain 

F  =  -16.33,  M  =  -3.77,  G  =  -0.95,  A'  =  -1.52,  L  =  -18.23,  N  =  -19.15. 
From  Eqs.  (37),  we  then  have  fCM”  the  effective  masses  in  the  <100> 
direction, 

mg  =  0.056mo  (0.058mo)(^o^ 
mj^P  =  0.36mo  (0.35mo)^^’^ 

m|j^  =  0.08mo  (0.09mo)^^’^ 

mgQ=0.16mo. 

The  values  in  parentheses  were  taken  from  the  references  indicated  and 
were  obtained  by  other  methods.  The  close  agreement  adds  validity  to  our 
model. 


The  primary  sample  (11316)  selected  for  presentation  in  Section  VII 
was  a  strained  layer  superlattice  composed  of  60  periods  of 

Inj^Ga^.^As/GaAs.  grown  on  an  lnoj6ao9As  buffer  layer  atop  a  <100> 


p-type  GaAs  wafer.  Well  and  barrier  thicknesses  were  120  angstrom  each, 
the  buffer  layer  was  Ijim  thick,  and  the  GaAs  barrier  layers  were  doped 
with  Si  over  20  angstrom  at  their  centers.  Strain  in  this  sample  was  0.7%. 


Carrier  concentration  was  9.6*10’*  cm'^  and  mobility  was  23,000  cmVv 
sec  at  4 
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VII.  DATA  AND  RESULTS 


A.  Hydrostatic  Pressure  and  flagnetlc  Field  Effects 

As  mentioned  in  the  discussion  of  samples  in  Section  VI.  we  selected 

Sample  t13l6.  an  Ir^  2^  strained  layer  superlattice,  to  show  the 

effects  of  hydrostatic  pressure  and  also  the  appearance  of  the  Landau 
levels  with  increasing  magnetic  field.  The  strain  in  this  particular  sample 
(0.7X).  due  to  the  lattice  mismatch,  lifts  the  valence  band  degeneracy  and 
presents  us  with  the  simple-band  situation  covered  by  the  calculations  at 
the  beginning  of  Section  V.  We  can  therefore  anticipate  the  energy  values 
of  the  Landau  levels  from  the  relation  (cgs) 

Efj  »  Eq  ♦  (n+l/2)ft(i>g  =  Eq  ♦  (n*I/2)licB/jJC, 
where  is  the  reduced  effective  mass. 

Since  the  pressure  coefficients  for  GaAs  and  other  Group  Ill-V  materials  at 
r  are  ail  approximately  0.012  *  0.002  eV/kbar,^’^^  we  can  expect  the  shift 

of  Eg  with  pressure  to  be  described  by 

Eg(P)  «  Eg(0)  ♦  0.012P. 
where  E  is  in  ev  and  P  in  kbar. 

From  the  combination  of  pressure  and  magnetic  field  data,  we  can  then 


determine  the  value  of  the  energy  gap  of  the  material  and  its  dependence  on 
the  external  parameters,  as  well  as  the  reduced  effective  mass.  )i.  and 
ApCP).  Figure  It  and  12  show  the  luminescence  spectra  at  a  pressure  of  1 
bar.  B=1  and  60  kC  respectively.  Although  data  were  taken  at  many  other 


values  of  the  magnetic  field,  these  examples  are  representative  of  the 
spectral  structure  observed  and  clearly  show  the  Landau  levels  at  the 
higher  field  value.  These  experimental  runs  were  then  repeated  with 
I  different  amounts  of  pressure  applied  to  the  sample.  Figures  13  and  14 

show  the  spectra  at  P  »  4  kbar  and  6  «  1  and  60  k6,  repectively.  In  these 
I  figures  the  logarithmic  intensity  has  been  added  (broader  curve,  right-hand 

scale)  since  the  Landau  levels  are  easier  to  discern  in  this  manner. 

Comparing  the  linear  spectra  at  P  =  I  bar  and  P  «  4  kbar.  and  noticing 
the  identical  shapes,  we  have  good  evidence  that  the  sample  was.  in  fact, 
subjected  to  almost  pure  hydrostatic  pressure.  Spectra  obtained  at  B»i  k6 

I 

I 

'  for  various  pressures  are.  except  for  the  peak  shift,  indistinguishable.  This 

is  an  excellent  validation  of  the  pressure  generating  technique.^^'^ 

An  indication  of  the  shift  of  the  luminescence  peak  energy  with 
pressure  is  shown  In  Fig.  15.  with  spectra  taken  at  P  =  1  bar.  •^2  kbar. 
^4  kbar  (left  to  right),  and  B  «  60  kC.  The  small  variation  in  these  spectral 


curves  resulted  from  the  fact  that  these  curves  are  from  different 
ruR'^ies  with  varying  power  levels  and  monochromator  slit  settings,  and 
do  not  indicate  non-hydrostatic  pressure  effects.  A  least  squares  fit  to 
the  peaks  as  a  function  of  pressure  gives  an  indication  of  the  pressure 
dependence  of  the  band  gap.  AEg  «  0.0122  eV/kbar.  This  is  shown  in  Fig. 
16.  We  can  thus  write 

Eg(P)«Eg(0)*  0.0I22P. 

where  Eg(0)  is  the  value  of  the  energy  gap  in  eV  and  P  is  in  kbar.  The 

measured  pressure  dependence  thus  fails  almost  into  the  middle  of  the 
expected  range. 

The  positions  of  the  Landau  level  energies  as  a  function  of  magnetic 
field  result  In  the  classic  fan-diagram’  shown  In  Fig.  17,  for  T  =  4  K  and 
pressure  of  1  bar.  The  numbers  associated  with  the  different  rays  of  the 
fan  indicate  the  particular  level  (nK>.1.2....).  Similar  results  were  obtained 
for  other  values  of  pressure  with  the  appropriate  shift  of  they -axis 

intercept.  Information  obtainable  from  these  plots  includes  Eg  and  the 

reduced  effective  mass.  A  good  approximation  of  the  energy  gap  of 
the  material  is  obtained  from  the  intersection  of  the  rays  ( n>2 )  with  the 
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y-axis  at  zero  magnetic  field  The  reduced  effective  mass  was  obtained 
from  the  slopes  of  the  lines,  where  we  considered 

I 

( 

I 

»  Eo*  (n*l/2)ft6>j.  =  (n+l/2)heB/cp. 

I 

Another  way  to  look  at  the  data  shown  in  Fig.17  is  to  plot  the  Landau 

I 

level  energies  as  a  function  of  ttw  normalized  magnetic  field.  (n«^l/2)B. 

( 

This  places  all  the  experimental  points  along  one  line  and  thus  provides  a  ! 

simpler  tool  for  analysis.  This  was  done  in  Fig.  18  for  P  =  Ibar. 

I 

The  intercepts  in  Figs.  17  and  18  show  Eg  %  1.318  eV.  As  a  first  order  | 

approximation  we  have,  for  the  direct  gap  in  lr\)  2^30  qAs, 

I 

Eg  «'l.52  -  I.60IX  ♦  0.536x2,  (68) 

or  Eg  1.2212  eV. 

To  this  value  we  must  add  the  changes  in  levels  due  to  the  wells  in  the 
conduction  and  valence  bands.  Using  the  effective  masses  from  the  k'p 
method  and,  for  this  approximatioa  assuming  infinite  wells, 

E,  =  (h27r2/2m"L/). 

where  m**  is  the  effective  mass  and  the  width  of  the  well, 
we  have. 


E,  c  *  0.0441  eV.  »  0.0069  eV.  E,  »  0.0309  eV. 


Since  we  are  dealing  with  a  conduction  band  to  light-hole  band 

transition,  we  obtain  as  an  estimate  of  the  energy  gap,  Eg  «  1.296  eV. 

This  calculated  value  differs  somewhat  from  the  experimental  result  (by 
‘^20mcv/),  but  is  within  the  range  of  variation  observed  across  samples. 

Thus,  if  at  the  point  of  investigation  is  100  instead  of  the  reported  120 

angstrom,  the  calculated  energy  gap  would  be  Eg  «  1.324  eV,  much  closer  to 
the  experimental  value. 

The  least-squares  value  for  the  slope  in  Fig.  18  is  p-'  =  20.7  ♦  0.2  (in 
terms  of  mo).  The  electron  effective  mass,  from  a  conduction  band  to 

acceptor  level  transition,  was  found  to  be  mg  =  0.069mo.^^5^  Using  this 

value  along  with  the  experimental  reduced  effective  mass,  we  can 
determine  the  effective  hole  mass  involved  in  this  transitioa  From 

l/)i  =  l/rr^  ♦  l/m^, 

we  have  m^  =  O.I6mo. 

Although  this  mass  is  larger  than  the  light  hole  mass  calculated  from  the 


k’p  method,  it  is  in  excellent  agreement  with  the  value  obtained  from 


oscillatory  magnetoresistance  measurements,  m^^  =  0.16mo.  Similar 


magneto-resistance  measurements  on  many  p-t|^e  samples  gave  results  for 
m^  ranging  from  0.13  to  0.17mo.(^s)  A  possible  explanation  for  this 

difference  between  the  calculated  and  measured  value  may  be  the  change  in 
the  band  structure  of  the  superlattice  versus  the  bulk  material.^^^^ 

In  Fig.  19  we  show  the  pressure  dependence  of  the  Landau  level 
energies  plotted  against  the  normalized  magnetic  field.  The  pressures  are. 
from  top  to  bottom:  3.83  kbar.  3.013  kbar.  2.005  kbar.  and  1  bar.  The 
reduced  effective  masses  may  be  obtained  from  the  slopes  of  the  fitted 
lines.  If  we  consider 

(mi/2)lie/M  =  AE/AB, 
and  write  >i  =  xmo,  then  we  have 

(n*l/2)l.l6xlO'Vx  =  AE/AB  (meV/k6). 
and  X  =  l,IB»10"  Vl/2)/£iE//:iB. 

For  the  pressures  listed,  we  obtain,  in  the  same  order. 

'-19.310.3, 19.710.2.  20.4i0.3.  and  20.7i0.Z 
Taking  the  two  extremes,  P  =  I  bar  and  P  =  3.83  kbar.  and  plotting  them  as 
in  Fig.  20.  we  see  the  pressure  dependence  of  the  reduced  effective  mass. 


<n  ♦  1/2)  B  <kG) 

Fig.  19.  Landau  level  energies  as  function  of  normalized  magnetic  field. 
Sample  n316.  at  T  =  4  ^K.  P  =  3.830  kbar,  3.013  kbar,  2.005  kbar,  and  1  bar 
(top  to  bottom). 
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Landau  level  energies  as  function  of  normalized  magnetic  field, 
316,  at  T  =  4  Ok,  P  =  1  bar  (upper)  and  P  =  3.83  kbar  (lower). 


Plotting  the  percentage  change  of  the  data  points  of  effective  mass  versus 
pressure,  as  in  Fig.  21.  yields  a  rough  estimate  for  the  pressure  effect  on 
the  reduced  effective  mass.  The  experimental  value  of  1.2  X/kbar  is 
comparable  to  the  calculated  change  of  0.8X/ld>ar(^)  For  this  calculation 
we  used  the  mass  terms  as  in  Eqs.  (37)  and  the  pressure  dependent 
constants  as  developed  in  Section  IV,  along  with  the  interpolated  values 
for  the  energy  levels.  We  obtained  the  conduction  band  effective  mass 
from  the  the  expression 

1  *  A'(ErEo)/(ErEo*AErAEo)  ♦  2P„2/3(Eo*AEo)  ♦  P„2/3(Eo*As*AEo). 
and  the  light  hole  effective  mass  from 

I  ♦  2FEo/3(Eo*AEo)  *  4GE2/3(E2*AE2)  ♦  ME,/3(Ei*AE,), 
where  the  constants,  including  P|^.  were  given  in  Section  III  and 

AEjsOEj/9P)jAP.  as  before.  Carrying  out  the  calculations  for  pressures 

of  1  bar  and  4  kbar,  we  found  a  percentage  change  for  the  reduced  effective 
mass  of  0.83*/kbar. 

We  discussed  earlier  that  in  undoped  GaA5/AI(Ga)As  structures 
transitions  are  exciton  dominated.  An  example  of  a  luminescence  spectrum 
showing  this  is  given  in  Fig.  22.  The  three  curves  shown  are  for  values  of 
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B  *  0.  27.530,  55.547  kG.  from  left  to  right.  This  sample,  T0082,  Is  a 
GaAs/AIAs  quantum  well  structure  with  60  Angstrom  wells  and  25 
Angstrom  barriers.  Data  were  taken  at  T  =  4  K,  P  =  I  bar.  In  addition  to  the 
quantum  well  transitions  shown,  this  sample  also  exhibited  luminescence 
peaks  near  1490  meV  (impurity)  and  near  1515  meV  (GaAs  layer).  Peak 
positions  for  all  three  transitions  are  plotted  against  magnetic  field  in 
Fig.  23,  with  a  curve  fitted  to  the  GaAs  transitions.  The  quadratic  nature 
of  the  shift  over  the  range  of  magnetic  field  investigated  is  clearly  evident. 
The  amount  of  the  shift  is,  however,  significantly  larger  than  indicated  by 

Eq.  (83).  Using  GaAs  effective  masses  as  calculated  we  obtain  Ry**  4.2 

meV  arxl  f1(Oo  ^  9.2  meV,  giving, 

At  (5/8Xl/2)^(fta)o)2/Ry"  »  0.8  meV, 


for  B  =  60  kG.  This  compares  to  an  observed  shift  of  2  meV.  The 
discrepancy  is  due  to  the  fact  that  Eq.  (83)  has  been  developed  for  a  fully 
confined  2-dimensional  excitoa  Exciton  energy  shift  due  to  a  magnetic 

field  however  depends  on  the  exciton  diameter  which  increases  with 

Our  experimental  results  agree  closely  with  those  obtained  by  others.^^^^ 


D.  Tmipcraturs  Effscts 


Since  the  primary  purpose  of  our  experiment  concerned  low 
temperature  (<4IC)  pressure  and  magnetic  field  effects,  temperature 
effects,  per  se.  were  not  investigated  in  detail.  However,  for 
completeness,  we  conducted  several  experimental  runs,  with  temperature 
as  the  variable,  to  determine  the  temperature  induced  shift  of  the 
luminescence  peaks.  This  was  accomplished  with  a  slightly  different 
experimental  configuration  than  described  in  Section  VI.  For  these 
measurements,  the  sample  was  mounted  in  a  vacuum  on  the  cold  head  of  a 
CTI-Cryogenics  closed  cycle  helium  refrigerator  and  illuminated  by  the 
laser  beam  through  a  window  in  the  Dewar.  The  refrigerator  was  capable  of 
providing  sample  temperatures  in  the  range  of  10  -  300K.  As  in  the  main 
experiment,  the  fluorescence  signal  was  directed  to  the  monochromator 
and  recorded.  The  results  for  a  selected  sample  (T00066).  a  GaAs/AIAs 
multiple  quantum  well  structure  with  ia^^r  thicknesses  of  40  and  25 
angstrom,  are  shown  in  Fig.  24  for  T  =  0.  120,  180,  240  K,  from  right  to 
left.  The  sharp  peaks  near  1515  and  1530  meV  are  grating  ghosts. 
Figure  24  clearly  shows  the  shift  to  lower  energies  with  increasing 
temperature,  as  well  as  the  growth  of  the  high  energy  tail  as  higher  levels 


are  thermally  populated. 


Applying  the  formula  for  the  gap  energy  as  a  function  of  temperature. 

E(T)«Eg(0)-ocTV(T*n  (71) 

with  the  6aAs  values  for  «  and  jl,  (  ot  =  8.871  *  10"^  eV/K,  ^  =  572  K),  and 
estimating  the  value  of  Eg(0)  from  Fig.  24  by  extending  the  straight  portion 

of  the  low  energy  slope  of  the  right-hand  curve  to  the  x-axis  intercept 
(*^1.572  eV),  we  obtain,  for  example,  for  240  K, 

E(240)  a  1.572  -  0.063  *  1.509  eV. 

This  is  the  same  value  obtained  when  the  low  energy  straight-line  portion 
of  the  left-hand  curve  is  extended  to  the  x-axis. 


I 


VIK.  CONaOSION 

The  experimental  investigation  of  different  semiconductor  structures, 
primarily  strained- layer  superlattices,  yielded  Interesting  results  which 
warrant  further  Investigatloa  Thus,  although  the  measured  pressure 

coefficients  for  Eg  fall  within  the  range  of  experimentally  determined 

values  for  the  appropriate  Group  fll-V  materials,  we  measured  consistently 
different  values  for  our  n-  and  p-type  samples.  Since  the  strain  in  these 
samples  has  different  values,  the  possible  effect  of  strain  on  the  pressure 
coefficient  should  be  investigated. 

A  similar  difference  was  determined  to  exist  in  the  pressure 
coefficients  for  the  reduced  effective  masses  of  the  n-  and  p-type  samples. 
Additionally,  our  measurements  indicate  significant  variance  of  the 
hole-mass  pressure  dependence  from  the  value  based  on  k  *  p  calculation 
for  bulk  material.  The  measured  dependence  is  considerably  higher  than 
anticipated.  A  possible  explanation  may  be  pressure- induced  valence  band 
mixing  in  the  superlattice.  This  aspect  should  be  investigated  theoretically 
and  through  further  experiments  with  specifically  tailored  samples  to 
determine  a  trend. 


The  developed  experimental  method  is  readily  adaptable  to 
measurements  involving  transmissloa  absorption  or  excitation 
spectroscopy,  and  experiments  should  be  conducted  to  Investigate  these 
topics  for  specific  samples  under  the  influence  of  pressure  and/or 


magnetic  field. 
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Appendix  A 
The  Empty  Lattice 

The  empty  lattice,  also  known  as  the  free  electron  model,  considers 
the  crystal  potential  to  be  equal  to  zero  but  maintains  symmetry 
properties.  IIl-V  compounds,  e.g.  GaAs,  form  in  the  zincblende  structure 
which  consists  of  two  interpenetrating  face-centered  cubic  lattlces.^^®*®) 
An  example  of  this  structure  is  shown  in  Fig.  The  First  Brillouin 

zone  of  a  face-centered  cubic  lattice  consists  of  a  truncated  octahedron 
and  is  illustrated  In  Fig.  A-2.  along  with  the  symmetry  points  and  axes.^®^ 
We  determine  the  eigenvalues  at  the  symmetry  points  r,  X,  L  of  the  reduced 
zone  for  the  face-centered  cubic  structure, 

^  =  (fi2/2mXk  ♦  6)2,  (A-l) 

where  6  *  reciprocal  lattice  vector.  The  primitive  translation 

vectors  joining  the  corner  of  the  cube  to  the  midpoints  of  the  adjacent 
faces  are 

Ti  =  (a/2)  (Oin  T2  =  (a/2)  (101),  Tj  =  (a/2)  (110), 
with  a  1  lattice  constant.  The  primitive  reciprocal  lattice  vectors  are 


obtained  from  the  equations  hj  =  2tt  (tj  x  ly)  /  (t,  •  12  x  *13). 


(A-2) 


where  i.  j,  k  are  cyclic  permutations  of  1. 2. 3.  We  obtain 

h|  =  (27r/a)  (111).  h2  =  (2ir/a)  (ill),  hs  =  (2Tr/a)  (ill). 

We  limit  ourselves  to  the  first  four  reciprocal  lattice  vectors  and  their 

variations  using  hj.  We  have. 

6o=  (2Tr/a)(000) 

G,=  (2Tr/a)(111) 

02=  (2Tr/a)(200) 

63  =  (2iT/a)  (220). 

The  variations  of  6j.  i  =0. 1,  2.  3.  are  the  following: 

60 :  (000) 

6, :  (111).  (TO),  (TO).  TOD.  (Hi).  TOD.  (iD),  (11D 

6; :  (200),  (020).  (002).  ODO).  (02D).  (0(2) 

6]  :  (220).  (202)  (022)  (&).  (22).  (089.  (22)).  (80).  202). 

(202).  0122).  (08). 

Using  these  values  of  6j  in  Equation  (A-l).  we  obtain,  at  the  point 


r  [k  =(2iT/a)  (000)1. 


En(r)  =  (Ii2/2m)  6)2  =  (li2/2m)(2TT/a)2  (000)2  I  =  0 

=  (li2/2mX2TT/a)2(111)2.  i  =  i 
-  (*2/2m)(27r/a)2  (200)2.  i  =  2 
=  (<»2/2m)(2TT/a)2  (220)2  j  =  3. 


At  the  point  X,  where  k  =  (27r/a)(10fl). 


Ejj  (X)  =  (Ii2/2m)  (2tr/a)2I(100)  ♦  O^R 


Substituting  for  6j, 


E|j(X.6o)  =  (fi2/2m)(2Tr/a)2(100)2 

E|((K6,)  =  (liV2mK2ir/aFC2llK  (OW.  CW,  (0H)*  (Oil)*, 


(011)2,  (2^,)2.  ciT) 


E|t(X,62)  =  (li2/2mK2lt/a)2(J00)2,  (120)2,  (io2)2,  (1)0)2,  (ia)2 


(1082 


Ek(X,63)  =  («i2/2mK2it/a)2CK0)2,  (302)2,  (122)2,  (a))2 


(12)2  (581)2,  (120)2.  (3,®2,  (1)2)2,  (12»2,  (ia)2 


At  the  point  L.  where  ki  =  (27r/a)(J  i  J),  we  have. 
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Eij(L.6o)  =  (l>V2m)(27T/a)2 ((iH)  ♦  6jl  or. 

Ek(L.Go)  -  (JH)2 

Ek(L.6,)  -  (3/2  3/2  3/2)2  (^2  V2 1/2)2  (3/2  V2  V2)2  (V2  3/2  V2)2 
(1/2 1/23/2)2.  (V2  3/2  3/2)2  (3/2  V2 1/2)2.  (3/2  3/2  1/2)2 
E|j(L.62)  -  (5/2  1/2  1/2)2  (1/2  5/2  1/2)2  (1/2  1/2  5/2)2 
(S/2  1/2  1/2)2.  (1/2  ^2  1/2)2,  (1/2  1/2  ^2)2 
£^(1.63)  -  (5/2  5/2  1/2)2  (5/2  1/2  5/2)2.  (1/2  5/2  5/2)2. 

(S/2  V2  1/2)2.  (^2  1/2  5/2)2.  (1/2  S(2  5/2)2. 

(5/2  ^2  1/2)2.  (^2  5/2  1/2)2.  (5/2  1/2  5/2)2, 

(^2  1/2  5/2)2.  (1/2  5/2  ^2)2.  (1/2  !^2  5/2)2. 

We  now  have  the  following  eigenvalues  with  degeneracies,  g.  as  indicated: 


r.g 

x.g 

ug 

(000)2  1 

(100)2.  2 

(1/2  1/2  1/2)2,  2 

(111)2.8 

(011)2,  4 

(3/2  1/2  1/2)2  6 

(200)2,  6 

(120)2.  8 

(1/2  3/2  3/2)2,  6 

(220)2,  12 

(211)2.4 

(5/21/21/2)2  &(3/2: 

(300)2&(122)2.5 

(5/2  3/2  1/2)2,  6 

v.v.v.v 


Appendix  B 


Fermi's  Golden  Rule  (so-se) 

We  consider  a  system  described  by  a  Hamiltonian 

which  satisfies 

i*(d/dt)|t>  =  H|t>  (B-1) 

where  Hq  is  time  independent,  but  H‘  may  be  time  dependent.  Then,  using  the 
interaction  picture,  we  define 

H'|(t)  a  (6-2) 

Substituting  Eqs.  (B-2)  into  Eg.  (B-l),  we  obtain 

ifi(d/dt)|t>i=H’,|t>|.  (B-3) 

Writing 

|t>|  =  |to>|*(HirJ,  *HW|f>,af. 
and  substituting  in  the  integrand  for  |t'>j^  then 


6 


These  substitutions  are  continued  until  we  arrive  at  the  expression 


|t>|  =U|(t,to)|to>  ]. 


(B-4) 


where  the  following  definitions  hold: 


» I 


(B-5) 


U|(t.to).  upon  using  Eq.  (B-2).  becomes 

UjCUoMl  Jt/-  e’^ot’/liH*(t*)e‘*^o(t‘-t’*)/l»H'(t‘')  * 

*  expl-iHo(t^J/fi)dt’dt"...dt  (B-6) 

Using  Eqs.  (B-4)  and  (B-6),  along  with  the  Definitions  (B-5),  we  can 
then  write  the  transition  probability  for  going  from  a  stationary  state  |to> 
1  I  i>  to  a  final  state  |  f>  after  a  time  t  as 


{B-7) 


Since  ( t,  to )  ■  I  results  in  <  t  |  i  >  =  0  (or  different  states,  the 


m 

m 


-INTI'S; 
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summation  starts  with  n  =  1,  and  for  first  order  yields 


P,f=|<f|Ul%to)|i>|2 


or,  Pjf  'ei(pli(Ef-E|)tV«l<f|H'(t')|l>dti  (B-8) 

Higher  order  contributions  to  the  transition  probability  may  easily  be 
developed.  Since  H'  is  usually  considered  to  be  a  small  perturbation,  only 
the  first  non-vanishing  term  of  Eg.  (B-7)  is  kept. 


If  we  define  the  operator  S|  as 


S, .  t.-«,U|(tto)  ■  SnSiH 


we  can  rewrite  Eg.  (B-7)  as 


Pir=l<f|s,|i>|2=  (B-9) 


This  allows  us  to  write,  for  the  first  order  contribution, 


P|f  =  I  I H'  I  i>ei(p(l(E  ,-E|)t/«l<lt  1 2.  (B-IO) 

Now  considering  the  Hamiltonian  to  be  of  the  form 
H*  =  K*e'‘“^  ♦  H’*e•“^ 

where  K'  and  its  hermitlan  adjoint  are  specified  to  be  independent  of  time, 
the  first  order  matrix  element  becomes 


<f  1  s'  1  i>  =  1 M'  1  i>  WtKE  f-Ej 


1 H'  I  i>  J.^em>(l(E,-Ej-*(i>)t/*l<lt, 
or.  by  using  a  definition  of  the  5-function. 

<f  I I  i>  =  -2Trl5(Ef  -E|-Ru))<f  1 K*  |  i> 


-2iTi5(E^-Ej*<ia>)<f  |  H‘  |  i>.  (B-11) 

where  the  first  term  on  the  right-hand  side  corresponds  to  absorption.  i.e. 
E(-Ej  s  lid),  and  the  second  term  to  emission.  Substituting  the  result  for 
absorption  into  Eg.  (B-10),  we  get 

Pjf  =  47T282(Ef-Ej-fia>)  I  <f  I H’  1 1>  1 2.  (B-12) 


where  52(Ef-Ej-lid>)  =  5(EpEj-fi<oX27tlir’ J_^**expIi(E  pE|-fi<D)t/lildt 

V 

I 


=S(E,-E|-*a)Xt/2W*)|.„”. 

Substituting  back  into  Eq.  (B-12)  and  defining  the  transition  probability  per 
unit  time  as  Wj^  >  (Pjf /t)  | 


we  arrive  at  a  form  of  the  "golden  rule",  first  order. 


Wjf  =  (27T/ft)8(Ef~Ej-fi4«))  I  <f  1  H'  I  i>  1 2  (B-13) 

In  order  to  determine  the  direct  absorption  coefficient,  we  write  the 
Hamiltonian  for  a  charged  particle  in  an  electromagnetic  field  as 

^total  ■  ^particle  *  ^radiation  *  ^  (B-14) 

where  H'  =  (e/2mc)(p  fl  ♦  fl  p)  ♦  (eV2mc2)A2  *  (e/mc)(fl  p). 

The  interaction  Hamiltonian  was  simplified  as  shown  by  neglecting  the 
squared  vector  potential  term  for  first  order,  and  using  the  fact  that 
(fl  p-p*fl)  can  be  shown  to  equal  zero.  Substituting  into  Eq.  (B-13),  we  then 
have 

Wj|  =  (2TT/li)(e2/m2c2)8(EpEj-liCi>)  |  <f  |  fl-p  |  i>  1 2.  (B-15) 

If  the  vector  potential  for  a  wave  of  frequency  a>  is  written  as 

fl(r,t)  =  ♦  c.c. , 

where  c  represents  the  polarization  direction,  and  we  recall  that  H'  = 
j^-g-iwt  then  we  have,  for  the  case  of  absorption. 

Wi(  =  (27T/fi)(eAo/mc)28(Ef-Ej-fia))  |  <f  1  e^'^'^ep  |  i>  1 2  (B-16) 

Using  the  dipole  approximation  on  the  matrix  element,  i.e.  expanding  the 
exponential  term  and  keeping  only  the  first  term  of  the  expansion,  we  have 


Wjf  « (27r/fi)(eAo/mc)28(^-Ej-fia))  |  <f  |  c  p  |  i>  p,  (B-I7) 


for  vertical  transitions  (kj=kf). 

Equation  (B-17)  can  be  used  to  determine  the  number  of  transitions  per 
unit  time  per  unit  volume  by  summing  over  the  possible  states. 

(d/dt)N||  ((i))«!(2iT/h)(eAo/mc)2 

«  Jg2*(2/8Tr3)  I  <(  |c.p  1 1>  I  28(e  f-E| -AO)),  (B-18) 

Where  the  summation  is  over  the  valence  and  conduction  bands,  and  the 
integral  Is  over  the  Brillouin  zone. 

Defining  the  absorption  coefficient.  wpCw),  as  the  quotient  of  energy 

absorbed  per  unit  time  per  unit  volume  and  energy  flux  (energy 
densityxvelocity).  then 

oc((ij)  =  ((d/dONjfCoi^fiwl/UxCc/n).  (B-19) 

where,  since  we  are  only  concerned  with  direct  absorption,  the  subscript  D 
will  henceforth  be  omitted,  and  the  energy  density  of  radiation.  U.  is  given 

t>y 


U  =  (n2/27Tc2)a>2Ao^. 


n=(real)lndex  of  refraction.  Combining  Eqs.  (B-18)  and  (B-19), 


ot((i>M47T2eVncm2a))2^^^Jg2^(2/8Tr3)  |  K'y  ^  I  *S(Ef-Ej-fio)).  (B-20) 

where,  for  convenience,  we  let  represent  the  matrix  element 
<f|c-p|i> 

Equation  (B-20)  shows  that  A((j>)  depends  on  to"’  and  an  integral. 
Removing  the  matrix  element  from  the  integral,  we  are  left  with  an 
expression  for  the  joint  density  of  states  for  a  pair  of  bands. 

Jc,v  -  Jg,(lt(2/8lT5)S(Ef-E|-ftO)).  (B-21) 


Using  Eq.  (B-21)  in  (B-20).  along  with  Eqs.  (37).  the  absorption  coefficient 
can  be  calculated. 
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